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BIOLOGICAL EPIMORPHISM, ADEQUATE DESIGN, 
AND THE PROBLEM OF REGENERATION 


N. RASHEVSKY 
COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


It is pointed out that the property of an organism of regenerating cer- 
tain lost parts is a subproperty of the general property of reproduction. 
The principle of adequate design imposes an upper limit upon the total 
number of subproperties which an organism can possess. Therefore with 
increasing complexity of an organism in some direction, which by defini- 
tion is the same as an increase of the total number of subproperties, the 
number of subproperties needed for regeneration must decrease. 


In a previous paper (Rashevsky, 1955) we attempted to explain 
from the point of view of the principle of biotopological mapping 
(Rashevsky, 1954) the well known fact that with the increase of the 
complexity and degree of differentiation of an organism, at least in 
animals, the ability to regenerate lost organs and parts decreases. 
The explanation based on the theory of point bases of graphs was 
not adequate. Its weakness becomes even more apparent when from 
a graph-theoretical approach we pass to a purely set-theoretical 
one (Rashevsky, 1958b, 1959, 1960a). A purely relational approach 
to the problem seems to be in general inadequate, inasmuch as a 
combination of the principle of epimorphism with R. Rosen’s in- 
teresting work on the problem of ‘‘reestablishment’’ of lost com- 
ponents in a relational system (Rosen’s J, R - systems) indicates 
the following: 

A component M of an ‘il, ® - system of greater complexity, which 
is constructed from a simpler ‘Il, ® - system according to the princi- 
ple of biological epimorphism, is reestablishable if the component 
My of the simpler system, on which M maps, is reestablishable 
(Rashevsky, 1960a). This definitely contradicts such facts, for ex- 
ample, as that an extremity of a salamander is reestablishable 
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whereas the corresponding extremity of a mammal is not, even though 
in Rosen’s terminology an extremity is, strictly speaking, not a 
component. 

The purpose of this note is to show.that while the principle of 
biological epimorphism, taken alone, may not be sufficient to ex- 
plain the decrease of regenerating ability with increasing com- 
plexity of an organism, this principle when combined with that of 
adequate design (Rashevsky, 1960b) offers a natural explanation 
which appears to be consistent with known facts. 

Let Np be the number of basic biological properties in the pri- 
mordial organism, while No is the number of all biological proper- 
ties and subproperties in a higher organism. We have always 


No>Np. (1) 


A particular biological property P; of the primordial may be carried 
by some more or less complex molecular structure. Any subproperty 
Pig, of P; in a higher organism is carried by a number of cells of a 
particular type, or even by a number of cells of different type. 
Thus the subproperty of vision, which is included in the property 
of general sensitivity (Rashevsky, 1958a, 1960b) and which maps 
onto the latter, is carried in a higher organism by cells belonging 
to a number of different types, which constitute the eye. 

The principle of adequate design requires that the number of 
types of cells, as well as the number of cells in each type, which 
carry a given biological subproperty P;«, must exceed a certain 
minimum in order to operate adequately. Thus a single muscle cell 
would be of no use in a higher organism if other muscle cells were 
not present. In fact, if the organism is to perform adequately the 
function of locomotion, the total number of muscle cells must be at 
least equal to a given fraction 7 of all the cells. If therefore we 
denote by \, the total number of cells in an organism, we have 


N. > No. (2) 


The fraction n;. of cells of a higher organism which carry a given 
subproperty Pi is partly determined, according to the principle of 
adequate design, by the nature of the subproperty P;« and by the 
total number N, of the cells, that is by the mass of the organism. 
It is also determined by the other subproperties which the organism 
possesses, as well as by the nature of its activities. Thus the 
relative number of muscle cells in an animal which leads a sedentary 
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life and which is protected from enemies will be less than in an 
animal which has to run around to catch its prey and has to ac- 
tively fight off other animals. 


Since 
2s Vie = 1s (3) 
i, & 


and since for a given type of organism each 1; must be equal to a 
minimal value, therefore the number No of all subproperties is 
limited by the design of the organism. Thus the principle of ade- 
quate design introduces a limitation on all the possible combina- 
tions of subproperties in a higher organism, even though this limi- 
tation is not formulated in an explicit and specific manner. This 
conclusion throws interesting light on the problems discussed in a 
previous paper (Rashevsky, 1958a; 1960b, chap. xxix), namely as 
to why we do not find for example organisms that perceive X-rays 
directly as some perceive visible light. We find here a different 
manner in which the principle of adequate design does limit the 
number of possible combinations of P;,’s. That such a limitation 
was needed and that it is likely to come from the principle of ade- 
quate design was emphasized in loc. cit. 

The above discussions also lead us directly to the conclusion 
that in higher organisms the regenerative ability must be smaller 
than in lower organisms. 

Amongst the basic properties P; of the primordial is contained 
the property Pr of reproduction. Amongst the subproperties Pro, of 
Pr we find such properties as the reproduction Pry of the whole 
organism, as well as the reproduction of the type of cells Cia 
which carry the particular subproperty P;q. A particular type of 
such a reproduction, in other words a particular subproperty Pra, 
is the production of a given cell type from other cell types, such 
as mesenchyme cells. The regeneration of a type of cell, or of a 
whole organ composed of several cell-types, is contingent on the 
reproduction of this type of cell either from cells of the same type 
or from cells of a similar type. 

If we have an organism which performs a large number of com- 
plex functions, in other words, an organism with a very high N,, 
then, because of the limitations introduced by (3), it cannot per- 
form some other functions. If an organism loses the subproperty 
Pry of reproducing itself, it will not multiply and will cease to 
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exist. The only subproperties of Pr, which higher organisms can, 
so to say, afford to lose are subproperties of regeneration of cer- 
tain of its parts. 

To put the argument in a different form, we may say that if the 
organism has a number N6 of subproperties other than regeneration 
of its parts, then, since each subproperty requires a fraction 77, of 
the total number N, of cells, only the fraction 


nr=1-)) te (4) 
i, O 


is left available for regeneration. The more developed and complex 
an organism, the higher, by definition, its Nj, and therefore the 
greater the sum of all n{.,’s in (4), and the less n,. 

In the light of the above discussion we can understand the re- 
sults of the experiments of Marcus Singer (1954) in which the re- 
generation of an extremity in a frog, which does not regenerate 
normally, is achieved by switching to the extremity some nerve 
supply from other parts of the body. This, as the author points out, 
indicates the need for a higher concentration of some special sub- 
stance at the site of regeneration. The production of this sub- 
stance, whatever its nature, is a subproperty P,q« of the general 
property P, of producing some biologically active substances. If 
this particular subproperty P,« serves only the purpose of stimulat- 
ing regeneration, then it can and will be weakened with increasing 
complexity of the organism. If it is produced by nerve fibers, as 
the experiments of Singer indicate, then the principle of adequate 
design may require for proper regeneration more nerve fibers than 
are needed for general nervous coordination. This may be incom- 
patible with (8) if a greater nervous coordination also requires the 
development of other subproperties. 

In connection with the problem of limitation of possible combina- 
tions of Pio’s, we should point out an interesting possibility which 
follows from some of Rosen’s (1960) discussions on the quantum- 
mechanical approach to genetics. Rosen (p. 253) suggests a rather 
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iconoclastic view that in the experiments on transformation of bac- 
teria through introduction of foreign DNA, the transformation oc- 
curs not because that DNA carries the genetic information but be- 
cause the DNA acts as the ‘‘observer’’ of genetic material carried 
by something else. The lack of a particular ‘‘observer’’ rather 
than of a particular information carrier is what produces, according 
to this view, the absence of a given phenotype. 

If this is so, then it may well be conceivable that in every or- 
ganism there is contained genetic information which, properly ‘‘ob- 
served,’’ would give many more combinations of subproperties com- 
patible with the principle of epimorphism than are actually ob- 
served. The limitation of the ‘‘observing’’ apparatus may be the 
reason for the nonoccurrence of many of those otherwise possible 
combinations. 


This work was aided by a Grant RG-5181 of the United States 
Public Health Service to the University of Chicago. 
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A CONTRIBUTION TO THE THEORY OF EGOISTIC AND 
ALTRUISTIC INTERACTIONS 


N. RASHEVSKY 
COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


If two individuals work on the production of the same object of satis- 
faction and if their cooperative efforts result in an increased overall pro- 
ductivity, then the following statement holds: if each individual tries to 
maximize his own satisfaction, each individual will have less of the ob- 
ject of satisfaction than when each tries to maximize the sum of the satis- 
factions of both individuals. This theorem was demonstrated previously 
for the case of a rather special choice of the satisfaction function. It is 
now shown that it holds for a much wider class of such functions. In 
addition to previously studied cases of methods of sharing the object of 
satisfaction, several new cases are investigated. An attempt is made to 
give a mathematical formulation for a society in which individuals con- 
tribute according to their ability and share according to their needs. It is 
shown that in order for such a description to be meaningful, certain re- 
quirements must be imposed on the satisfaction function. Even with 
those requirements fulfilled, the possibility of such a society requires the 
fulfillment of an inequality between certain functionals. 


In several previous publications (Rashevsky, 1947a, 1947b, 1959) 
we investigated the interaction of two individuals who cooperate in 
the production of some ‘‘object of satisfaction.’’ We considered 
two basic possibilities: first, when each individual tries to maxi- 
mize his own satisfaction function regardless of what the other 
does (egoistic or individualistic behavior), and second, when both 
individuals try to maximize the sum of the satisfactions (altruistic 
or collectivistic behavior). We assumed that the satisfaction func- 
tion increases logarithmically with the object of satisfaction, 
that the work or effort necessary to produce the object of satisfac- 
tion creates a negative satisfaction proportional to the amount of 
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the work or effort, that the productivity of two cooperating indi- 
viduals is greater than the sum of productivities of each individual 
working separately, and that the total amount of the object of satis- 
faction produced by both cooperating individuals is divided between 
the two individuals either in proportion to their individual effec- 
tiveness or in a fixed ratio. Under those conditions we found that 
the altruistic behavior leads to the production of greater amounts of 
the object of satisfaction and to greater individual satisfaction. A 
number of interesting studies somewhat along the same lines have 
also been made by A. Rapoport (1947a, b, c). 

The purpose of this paper is to study a few other cases of this 
type of interaction. Before doing this we shall, however, generalize 
somewhat our previous results and show that they do not depend on 
the particular choice made previously for the form of the satisfac- 
tion function. 

I. We denote by z the amount of the object of satisfaction and by 
S(x) the satisfaction due to that object. We assume that for z > 0 


ads d785 dS dS 
(x) a Pie : (45) ; (35) 0, -,d) 


and such that dS/dz decreases faster than 1/2. 

As in the previous papers we denote by y the work necessary to 
produce the amount z of the object of satisfaction, and denoting by 
B a constant, we consider that the negative satisfaction due to y is 
equal to- By. Hence the total satisfaction Sp is equal to 


Sr =S(x) - By. (2) 


Also as before we assume that z is proportional to y. If @ is the 
factor of proportionality, then 


v= ay. (3) 
Hence from (2) and (3) we have 
Sr(y) = S(ay) - By. (4) 


Because of (1) Sy has one and only one maximum for a positive Yy- 
Denote by the subscripts 1 and 2 the symbols which correspond to 
two different individuals and consider the case where two indi- 
viduals cooperate in the production of the same object of satisfac- 
tion. Then, as before, we assume that the total amount X of the 
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object of satisfaction produced by the two individuals working to- 
gether is 


X=@iyit+@oy2+a3yiyo- (5) 


If the amount X is divided between the two individuals in prupor- 
tion to the effectiveness of their work, then the corresponding shares 
are 


Q@1Y1+ @2Y¥2+ @3yY1Y2 


v1, = RH tikes Vs Q1Y15 (6) 


@1Yi1+ @2Y2 + 3 Y1 Y2 


ras 21Yi+ @2Y2 eee ae 7) 


Equation (2) holds for each individual whether working sepa- 
rately or in cooperation with the other. In the first case because of 
(3) we obtain for each individual expression (4). In the second 
case we have 


Sr, = 81 (@1)-Biyi; 
Spo = S2(v2) - Baye. 


Substituting (6) and (7) into (8) we obtain S7,; and Sr2 as func- 
tions of y; and yo. 
For egoistic behavior the values y; and y2 are obtained from the 


(8) 


equations: 
fs] 0s 

Side Faye BERTH Tey (9) 

Oy 0 Y2 

Put 

as Sr 
Tt = Fi(yis y2)3— = = F2 (yt 2) (10) 

OY: J Y2 


and consider the curve F (y1,y2) = 0 
We have from (8): 


Fi(yisys)= Set = StS 8, a 


From (6) we see that for y2 = 0 
fe) v1 
= 2 ae (DAN, (12) 
C1, =41Y1; ay, 1 
Hence because of (1) and (11) the curve F; (y1, y2) = 0 intersects 
the y;-axis at a finite point yi = C which is determined by a1, Bi 
and any other parameters which may enter into S; (x 1). 
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Now consider the case y2 =». We have from (6): 


a1aoy1+41a3y17 | Ox, _ 41 Ga + 201 ag yt (13) 
Si a2 Oy, ao 


Hence for y2 = ~ the derivative dS;:/dz; becomes a function of 
y1 only. So does d21/dy1. On the other hand it is seen from (6) 
that dv,/dy; is always positive while according to (1) dS;/dzx 
decreases with x; faster than 1/z;. For very large values of y1, 
a, increases as 71, while d2;/dy1 remains constant. Hence the 
expression 


decreases faster than 1/y,; and therefore 


dS dep, 
dt, Oy, 

is negative for very large values of y;. For y; = 0 we find from 
(18) that x, =0, dxi/dy1 > 0. Hence, because of (1), the above 
expression is positive fory,; = 0. Because of (1), (dS:/dz1)(0z1/dy; ) 
decreases monotonically. Hence for y2 = ~ the equation F; (y1, y2) = 
0 has a finite positive root y; = y;- This means that the curve 
F';(y1,y2) = 0 has an asymptote parallel to the yo-axis. 

By a similar argument we find that the curve F2 (y1, y2) = 0 has 
an asymptote parallel to the y;-axis. Since the asymptotes inter- 
sect, therefore the curves F; (y1, y2) = 0 and F2(y1, y2) = 0 inter- 
sect in at least one point (y;*, yo*). 

Since nothing has been said about the monotonicity of those 
curves or of their derivatives it is conceivable that they will inter- 
sect at more than one point, as shown in Figure 1. The number of 
intersections is always odd. We shall return to this case later and 
consider now the case where the two curves have only one point of 
intersection as shown in Figure 2. 

As before we find from (6) and (7) 


On ajay rs) 2 a3 y.2 
Aes 1 7341 > 0; iit es eR 5>0. (14) 
Oyo (a1 y1 + ag y2) Dy: (a1 y1 + a2y2) 
Consider now the case of altruistic behavior. Here the values of 
y1 and yo are determined by 
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Fly, yo) 


Ye F,(y,,¥o) 


*yl) yon 
y, PA: 


FIGURE 1. 


0 (Sr; + Sra) ee 0 (Sr, + Sre) ~ 


. OF 
Oy; IY 
This may be written: 
OSr, aAS8re OSt,; ISre2 
+ = 0; Se nS 0. 
OY, Oy; dY2 dY2 
F, (y,,¥2) 


Fay), Yo) 


FIGURE 2. 
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But because of (1) and (14): 
OST, aST2 O22 


OS (17) 


OY, Ax » OY1 


Therefore d(Sr; + Sr2)/dy1 > 0 at all points of the curve 
OSr,/dy1 = 9. Since for a fixed finite yo, as yi — ~, % remains 
finite, while dz2/dy1 decreases to zero as seen from (14), there- 
fore OSp/dy, tends to zero as yi—> ~. Hence the locus of the 
maxima of Sr, + Sr with regard to y1 will be to the right of the 
line F; (y1,y2) = 0. (Figure 2). By a similar argument we see that 
the locus of maxima of Sy; + Sr, with respect to y2 will be above 
the line F'2(y1,y2) = 0. Hence the values y: and yo for altruistic 
behavior are larger than for egoistic; so will be the values x; and 
#2, just as we found in our previous papers in which a special form 
of S(z) was assumed. 

The restriction (dS/dz),29 = ~ in (1) is not necessary. However, 
if we do not make it and still require that Sy (y) in (4) have a maxi- 
mum for each individual, it would be necessary to add the re- 


x=0 


We shall now consider the case represented in Figure 1. For 
egoistic behavior we now have three solutions of the equation (9). 
It can, however, be seen that the solution which corresponds to the 
point O» is unstable. Consider a situation which corresponds to 
the point y,“?, y{Y in Figure 1. Since the curve F, = 0 represents 
the locus of maxima of Sr, with respect to y;, and since the first 
individual tries to maximize his satisfaction Sy, with respect to 
y1 for any y2, therefore the first individual will try to decrease y;. 
By the same token the second individual will try to increase yo. 
As a result the point y“¥, y{? will be shifted towards 0;. By a 
similar argument we see that a point (y*, y{) will be also shifted 
towards O,. A similar situation prevails near the point O3;. Thus 
O; and Og represent stable states. But as seen from Figure 1, any 
displacement from Og towards either 0; or 03 will increase. Thus 
Oz represents an unstable state. Therefore only either the state 
which corresponds to 0; or the one that corresponds to O3 will be 
possible. A change from egoistic to altruistic behavior results in 
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a shift of the curve F’; = 0 to the right and of the curve F2 = 0 up- 
ward. Such a displacement moves both 0, and Qe to the right and 
upwards. Thus, regardless of whether we consider the state which 
corresponds to 0; or the one which corresponds to 02, we see that 
altruistic behavior results in an increase of both y; and yo. If the 
parameters of the different functions involved are such that the 
shift of the curve F; (yi, y2) =0 to the right and the shift of the 
curve F2(yi1,y2)=90 upward are sufficiently large when passing 
from egoistic to altruistic behavior, then the situation shown in 
Figure 3 will occur. The points 0; and QO» disappear and only O03 
remains. If we start with the state O; for egoistic behavior and 
pass to altruistic, we have as before an increase of y; and yo. If 
we start with the state 0; for egoistic behavior we shall have an 
even more pronounced increase in y; and yg. Thus our results hold 
in this case also. 

The curves F; (yi, y2) = 0 and F2 (y1, y2) = 0 may intersect their 
asymptotes, as shown in Figure 4. More complicated cases may 
then arise. Perhaps it may be proved that conditions (1) are suf- 
ficient to insure the monotonicity of F; (y1, y2) = 0 and F2 (y1, y2) = 0. 
If so, then the situation discussed in the preceding paragraph can- 
not occur. Until, however, such a proof is given, we must con- 
sider the possibilities illustrated in Figures 1 and 4. 


Ye 


Os 


) % 


FIGURE 3. 


122 N. RASHEVSKY 


Yo 


y, 


FIGURE 4. 


Il. Now let us consider the case of division of the total amount 
of object of satisfaction in a fixed ratio (Rashevsky, 1947a, b; 1949, 
1951). Then we have: 


v = B(a1y1 + G2 Y2 + G3 Y1 Y2)3 (19) 
@ = (1- B)(aiy1 + a2 y2 + a3 y1Yy2)- (20) 


By an argument similar to the one used in the previous section 
we can show that the conclusion of the previous section holds 
again, thus generalizing our previous results. 

We may ask what determines the fraction 8B? In altruistic be- 
havior when both individuals try to maximize the sum Sr, + S72 of 
the satisfactions, they will agree to choose such a value of B which 
gives the largest value of Sp, + Srg. In other words, to the system 
of equations (15) there will be added the equation 


0(Sr, + Sto) 7 

dB i 

The system (15) and (21) contains three equations for the de- 
termination of y;, y2, and B. 


As an example we may take the particular form of Sy () used in 


loc. cit. We have in our present notations with A;, Ae, 41 and &% 
as constants: 


0. (21) 


Sr, = A; log a1 B21-Biyi; (22) 
St, = Ag log d2(1- B)w2-Biye. (23) 
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The equations 


O(Sry + Sto) ere O(Sry + Sto) ae 4 
Oy; OY¥2 ‘. 


? 


are in this case independent of 8 and determine the values of Y1 
and y2. Equation (21) now gives: 


Ce) 
ypu log B + Ag log (1- B)] = 0, (25) 
or 
A Ag 

= = 0, (26) 

Bias 
This expresses £ in terms of the parameters A; and Ag, namely: 

A 
ae 27 
B A, + Ao ( 


That this value of 8 corresponds to a maximum is seen from the 
fact that the second derivative of the expression in brackets in (25) 
is 

Ay Ag 


pete Ay 

How about the case of egoistic behavior? In this case it is to 
the advantage of the first individual to make £ as large as possi- 
ble, while it is to the advantage of the second individual to make it 
as small as possible. On what value of f will they compromise? 

Consider a fixed value of 8. For this particular value the opti- 
mal values y; and y2 are expressed in terms of 6 and of the other 
parameters which enter into the problem. Therefore the maximum 
values of Sr, and Sr. are functions of 8. We denote them by 
S7,(8) and St2(f). Denote by S;, and S;, the maximum values of 
Sr, and Sr, for the case where the two individuals work sepa- 
rately. Consider the differences 


581 = 87, (B) - 8,3 (28) 
582 = ST2(B) -Sr.- (29) 


If 8 is very small, then 5S; may be negative. In other words, in 
spite of the increased total amount of the produced object of satis- 
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faction as given by (5), the first individual in this case gives away 
so much of the product to the other, that he does actually lose as a 
result of cooperation. Excluding physical coercion on the part of 
the second individual, a situation which we shall not consider here, 
the first individual will not cooperate at all for such a low B. Simi- 
larly, too large a B will make 6S, <0 and the second individual 
will not cooperate. Thus the requirement 


5S; >0:— 68950 (30) 


limits the possible range of 8. If within that permissible range of 
B we have 5S, > &Sg, then if the second individual refuses to co- 
operate he will lose less than the first. Similarly if 6S; < 6Se, 
then if the first individual refuses to cooperate he will lose less 
than the second. Thus in the first case the second individual can 
force the first to a decrease in B by refusing to cooperate. In the 
second case the first individual can force the second to an increase 
in B. 

When 6S, = 6S2, both individuals stand to lose equally by re- 
fusal to cooperate. Hence the value of 8 on which a compromise 
will be reached is given by the equation 58; = 5S2. Thus 8 is de- 
termined by 


St,(B) - Sry = St2(B) ~ Spa. (31) 


Another possibility which may be considered is that both indi- 
viduals agree to divide the object of satisfaction in such a manner 
that the satisfaction Sr; and Sr2 of both individuals are equal. 
Since Sy; and S72 are functions of yi, y2, and #, therefore the re- 
quirement of the equality of Sp, and Srg leads to 


Sri (yi, y2,B) = Sra (yi; y2; B). (32) 


This equation determines 8 as a function of y; and ye for any 
values of the latter. Therefore Sp; and Srg now become functions 
of y; and yo only. From this point the individuals may either be- 
have egoistically or altruistically. In the first case the values of 
y1 and y are determined from (9); in the second case, from (15). 

III. In previous publications (Rashevsky 1947a, b; 1959) we gen- 
eralized some of the results obtained for one individual for any 
number of individuals. In that case we can speak of egoistic and 
altruistic societies, or of individualistic and collectivistic socie- 
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ties. A further generalization is made by considering intermediate 
cases, in which each individual maximizes a linear combination of 
the satisfactions of all the individuals. Egoistic behavior corre- 
sponds to the case in which the linear combination for each indi- 
vidual consists of his own satisfaction only. Altruistic behavior 
corresponds to the case in which for every individual the linear 
combination consists of the sum of satisfactions of all individuals 
(Rashevsky, 1947a,b; 1959). In the general case the behavior of 
society is characterized by a matrix (Rashevsky, 1959) which 
Anatol Rapoport appropriately calls the ‘‘social matrix.’ In an 
egoistic or individualistic society the social matrix is diagonal. 
In an altruistic or collectivistic society all elements of the matrix 
are equal to 1. Intermediate types of behavior are represented by 
more general types of matrices. 

We have called attention (Rashevsky, 1959) to the circumstance 
that the mathematical models of egoistic and altruistic societies 
show traits which bear some resemblance to socio-economic forms 
which are widely discussed both practically and theoretically. Al- 
though the actual situations are considerably more complicated, the 
abstract models may be very useful in clarifying ideas which are 
of great practical importance. We suggested beniatively, (Rashevsky, 
1959, p. 147) that we may conveniently define ‘‘economic individu- 
alism’’ or ‘‘free enterprise’’ as a socio-economic system which 
corresponds to an egoistic society in which sharing is made on the 
basis of efforts involved. ‘‘Socialism’’ may be defined as a system 
which corresponds to altruistic behavior, with sharing again made 
on the basis of efforts involved. To the extent that we are free to 
define any term the way we want, the above definition, which boils 
down to a mathematical description, is no worse than a number of 
much less precise definitions used by social and political sci- 
entists. It has the advantage of definiteness and its acceptance may 
well lead to a better mutual understanding. H. G. Wells (1949, p. 975) 
once wrote to the effect that every reasonable person is a socialist but 
nobody knows exactly what socialism is. While the factual correct- 
ness of the first part of the statement may be subject to some doubt, 
it seems fairly safe to subscribe wholeheartedly to the second part. 
Definitions are basically arbitrary. Yet useful definitions are to be 
preferred and a mathematical definition is likely to be more useful 
because of its precision and because of the possibility of drawing 


clearcut conclusions from it. 
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From the definition suggested above, we see at once that there 
is a continuous transition between ‘‘economic individualism’’ and 
extreme socialism. Inasmuch as the individual elements in the 
social matrix can vary continuously within a certain range, the 
total number of possible social matrices and hence of possible 
social structures is not only infinite but even not countable. 

It would be naive to assume that an actual socio-economic struc- 
ture is completely described by the social matrix. While given ac- 
tual structure, as we said above, corresponds to a definite social 
matrix, other characteristics are necessary to describe an actual 
society completely. This, however, does not reduce the usefulness 
of abstract mathematical formulations which deal with oversimpli- 
fied situations. 

From this point of view it is worth investigating what mathemati- 
cal definition may be appropriately given to a still utopian but 
widely discussed social structure, namely the Marxian communist 
society. It is usually described as a society in which every one 
contributes according to his ability and receives according to his 
need. While ‘‘ability’’ in this case may naturally be understood as 
efficiency or productivity of work, the term ‘‘need’’ has not been 
properly defined in any biological or socio-economic terms. Very 
tentatively we have suggested (Rashevsky, 1959, p. 147) that a 
communist society may be defined as a collectivistic society with 
sharing determined by the requirement of equality of everyone’s 
satisfactions. Such a society would be a generalization of the 
case of two individuals mentioned at the end of section II of this 
paper, and characterized by equation (32). This, however, does not 
seem to be an appropriate definition, although a society of this 
type may well be a legitimate object of mathematical study. 

When we feel a need for something our satisfaction increases as 
we obtain that something. As long as an increase of the amount of 
an object of satisfaction increases our satisfaction we feel some 
need for it. Thus the need may well be measured by the derivative 
dS/dz. This suggests the consideration of such a method of shar- 
ing in which, in the case of two individuals, dS;/dx; = dSo/dao. 
We again can introduce the fraction 8, determine for a fixed the 
equilibrium values y;* and y.*, find from (19) and (20) the corre- 
sponding values of x; and wa, and determine the values of dS;/dz 
and dS2/dx_ for those values of the variables 2; and x3*. The ex- 
pressions (dS;/da1),,=x,* and (dS2/dz2)x,=x,+ Contain B as param- 
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(=) aSe 
day ee dy %o™= xo" 63) 
now determines B. 


In such a case both individuals will have the same need for x, 
and it would hardly make any sense to talk in this case about a 


eter. The equation 


distribution according to need. The italicized words imply a dif- 
ference of needs. 

The only possibility left seems to be to consider such a distribu- 
tion in which the function S$; (x;) and S2 (x2) would have a maximum 
for each individual. But as long as we consider that S (7) increases 
monotonously with z this is impossible. Under those conditions 
the need of an individual is never satisfied. True enough, the total 
satisfaction Sr has a maximum. But this maximum is due to the 
negative satisfaction — By, caused by the efforts. This maximum 
does not reflect the actual satisfaction due to the possession or 
use of z, the object of satisfaction. 

We have pointed out, however, (Rashevsky, 1959) that in many 
instances, the satisfaction function S(z) has a maximum. A certain 
amount of the object of satisfaction is the most pleasant. Not 
enough of it reduces our satisfaction; too much of it also reduces 
it. In fact this situation may be more general than the case of a 
monotonically increasing satisfaction. 

If S(z) has a maximum for 2 = 2*, then it is natural to say that an 
individual needs an amount z* of the object of satisfaction. If an 
individual can obtain the amount 2* without any work or effort, say 
as a gift, then his need can be fully satisfied. If, however, he has 
to perform an amount of work y in order to obtain a, if again @ = ay, 
and if y results in a negative satisfaction - By, then, as we shall 
now show, the total satisfaction S; of the individual either does 
not have a maximum at all, or has a maximum for a value z = 2° 
such that 2” < 2*. 

To see this we first express the negative satisfaction in terms of 
z. It is equal to -(B/a)z. If (dS/dz),=9 < B/a, then unless S (z) 
has inflection points, so that dS/dz first increases, then decreases, 
the function Sr (zx) = S(x) - (B/a)z will have no maximum at all. If 
(dS/dz),= 9 > B/a, then S7(z) has a maximum. However, since for 
2 = 2* we have dS/dz = 0, therefore for 2 = x* we have aSt/dx = 
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dS/dz - B/a <0. This means that the maximum of Sr corresponds 
to such a value 2’ of x that x’ < a*. 

Hence, as long as the individual has to perform some amount of 
work, whether he works alone or in cooperation with another in- 
dividual, and as long as work in any amount is unpleasant, the in- 
dividual can never satisfy his need completely. 

Work, however, does not need to be unpleasant in any amount. 
Although there may be pathological individuals who refrain from any 
activity, yet in general a healthy individual will like and need to 
perform a certain amount of work. True enough, he will in general 
like some kind of work and dislike another kind. The existence of 
work of different kinds will have to be considered in a further de- 
velopment of the theory. At present we shall limit ourselves to the 
abstract oversimplified case, in which we can speak of “‘work in 
general.’’ This is no greater a limitation than the consideration of 
only one object of satisfaction. In effect we are considering one 
object of satisfaction and the one type of work which is involved in 
its production. 

If to some extent work is pleasant, then the satisfaction s (y) 
due to it is not always negative. It has a maximum for a value 
y=y*. Since z = ay, therefore we have s(y)=s(2/a). If the in- 
dividual performs the amount y* of work, he will obtain a maximum 
satisfaction from this work. Beyond the point y*, the satisfaction 
from the work decreases and may even become negative. 

Consider now the total satisfaction 


Sr(z) = S(x) +s (:). (34) 


Since for = 0 we have dS/dz 2 0, ds/dz > 0, while for suffi- 
ciently large values of z we have dS/dz < 0, ds/dz < 0, therefore 
at least for one « = 2** > 0 we have dS/dz + ds/dz = 0 and this 
corresponds to a maximum of Sy. Thus Sr has at least one maxi- 
mum but it may have two or more. 

If the optimal value y* of work which maximizes s is not enough 
to produce the optimal amount 2* of 7, so as to maximize S, then 


ay* < oF; omt/a>igts (35) 


d 1 
as) : (22) <0. (36) 
de x= x* a dy y=x*/a>y* 


Wence 
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Since (dS/dz),.,+ = 0, therefore it follows from (36) that 


es) dS ds 
Ae ane. = a aoe + Try East <a Oe (37) 


Hence the maximum of Sr lies to the left of e*. In other words, 
if «** denotes the value of x for which Sr has a maximum, we have 


ee Sees (38) 
Similarly we show that when (35) holds and if we denote by x’ the 


value of x for which s(y) = s (= has a maximum, we have 


z** > 2’, (39) 


Denote by Sp the maximal value of Sr, while by S* and s* we 
denote the maximal values of S(x) and of s(z/a). We have, because 
of (38) and (89): 


St = S(x**) + s(a**/a) < S* + 8*. (40) 
Only when 
ay* i r*, (41) 


that is when the optimal value of y* is just enough to produce z* 
do we have Sf = S* + s*, and this is the largest value which Sp 
can assume. 

If ay* > x*, we show in a similar manner as before that 2** > 2*, 
2** <2’. Therefore in that case equation (40) still holds. If we 
assume that the excess of x over x* produced in this case can be 
stored, and that only the amount z* will be used by the individual, 
then for ay* > z* a single individual can simultaneously maximize 
both his S and his s, and thus have his ‘‘need’’ for both the object 
of satisfaction and for work completely satisfied. Such a case is, 
however, much less likely to occur and presents less interest than 
the case represented by inequalities (35). In that case, an indi- 
vidual working alone cannot completely satisfy his need. 

If, however, two individuals cooperate, then the situation may 
change. We may now consider a more general expression than (5). 
We may put 

X =ary1 + aoy2 t+ u(y1s¥2)s (42) 


where the function u(y1,y2) is defined for y, > 0, y2 > 0, and sat- 
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isfies the following conditions 


Ou Ou 
‘ = 0) =0; (——»>.03) 2 Dae 
u(yi,y2) > 9; u (0, v2) vACRE, ) ’ ai dys ( ) 
We shall again denote with the indexes 1 and 2 correspondingly 
the different symbols which refer to the first and second individual. 
Now, even if a; y* < ai, and a2 y2 < #*, a situation which corre- 
sponds to (35), it may happen that 


Q1yi + agye+ulys,ys)=2r +e . (44) 
This requires that 
u(yr,y2) =a + ao -ayy1 —adeyo- (45) 


Under those conditions the first individual can use the amount 
z;* of z, the second the amount z,* of z, while performing corre- 
sponding amounts of work y;* and yg*. Thus each will work accord- 
ing to his best desire and receive according to his need. 

If 


Q1yi + agyo + ulyr, yo) > 2 t+ ae", (46) 


and if again the excess of X over wy + 2g" can be stored or given 
away, then under those conditions a distribution according to need 
is again possible. If, however, 


Ary + aoy2t+ulyrs yo) <r + 2o%, (47) 
then it is impossible to satisfy simultaneously the needs of both 
individuals. 


Thus the possibility of satisfying the needs of both individuals 
simultaneously requires that 


u(yt, yo") S ai* + @a° -—aiyt + aye’, (48) 


which means that the cooperative effort must yield a sufficiently 
large increase of the total production output. 

The ‘‘need’’ of an individual is expressed by x*; his desire to 
work by y*, his ‘‘ability’’ by a. If it happens that for one of the 
individuals, say the second, yo" = 0, ze* > 0, and if for the first in- 
dividual a, y;* = x\* + xo", then the first individual will maximize 
his Sp by putting in an amount y;" of work, producing an amount 
v;* +29" of z, and giving away the amount 23" to the second in- 
dividual. This is still possible if a; y* > 2;* + ao* and the excess 
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of x over 7;" + 2" can be stored. Thus each individual will have 
his needs completely satisfied but the second will be parasitizing 
on the first. 

We can generalize the situation to the case of No individuals. If 
No is very large, we consider the distribution function N (@*,cytageer)s 
so that N (2*,y*, a)dx*dy*da gives the number of individuals who 
have an z* between x* and x* + dx*, a y* between y* and y* + dy*, 
and an a between a and a + da. 


We have 
| i | N (x*, a, y*)dx*dy*da =N,. (49) 
0 Yo ~o 


The function 


Ni(a, y*) -| Nie*,a*, a) dz* (50) 

) 
gives us the number N, (a, y*)dady* of individuals with a given y* 
and a. The function 


Na (e*) = [ N, (2*, y*, a) dady* (51) 
0 re) 


gives us the number of individuals N» (2*) with an x* between zx* 
and #* + dz*. 

For each individual working alone we have z = ay. Therefore if 
all individuals were working separately and each put in his optimal 
amount y* of work, then the total amount of x produced by the N in- 
dividuals would be 


$a ah [ N, (a, y*) ay*dady* . (52) 
0 Jo 


When working in cooperation, the N individuals will produce an 
amount X > X’.. The excess X - X’ will depend on all the a’s and 
all the y*’s and thus be a functional F'[N; (a, y*)] of N1 (a, y*) just 
as X’ is a functional of that function. Thus we have 


X= ihe if N 1 (a, y*) ay* dady* + F (IN: (a, y*)]- (53) 
0 


0 
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The total amount of 2 necessary to satisfy everybody’s need of 


it is given because of (51) by 
Xq* -| 2*No (2*)de*. (54) 
f) 
If 


0) 0) 


0 


then a satisfaction of everybody’s need is possible. In a case 
where the inequality sign in (55) holds, everybody’s satisfaction is 
possible only if the excess of the total-2 can be stored, provided 
the process of storing does not reduce the satisfaction. In order 
that inequality (55) could be satisfied, a sufficiently large produc- 
tivity is necessary. If we tentatively discuss this highly abstract 
and oversimplified picture in real terms, we may say that inequality 
(55) means a high degree of industrialization both in regard to out- 
put and proper organization. 

In general there will be a small but finite number of individuals 
who have a finite or even possibly a large need z*, but who also do 
not like to work and whose y* will be practically zero. The total 
number of such individuals is obtained by considering the function 


Na (2*,y*) = | N (e*,a,y*)da (56) 
0 


which gives us the number N3 (x*, y*)d2*dy* of individuals whose 
z* lies between x* and 2* + dx*, and whose y* lies between y* and 
y* + dy*. The number of individuals with a very small y* < § and 
any given z* is 


Ng = 5 | Ng (x*, 0) dx*, (57) 
ft) 


and the fraction of such ‘‘parasitic’’ individuals is N4/No. 

The existence of such ‘‘parasitic’’ individuals may be considered 
as a serious disadvantage of the type of society considered here. 
If due to imitation effects (Rashevsky 1959) their number increases, 
this would make the type of society considered here unstable. How- 
ever, in any society which existed in the past or exists now, there 
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always has been and is a fraction of individuals who are ‘‘parasitic”’ 
in the sense that they cannot take care of themselves and must be 
supported by the rest of the society. Cripples, mental defectives, 
etc., belong to that class. It is more than likely that the normal 
functioning of the healthy human organism requires a certain amount 
of work, an amount which is pleasant. In that case the fraction 
N4/No may play the same role in our hypothetical society as the 
fraction of individuals that cannot support themselves plays in real 
societies. Moreover, the desire or love for work is likely to be to 
a large extent culturally conditioned, something that cannot be 
said about many purely physical requirements necessary for work in 
actual societies. All the above points offer interesting problems 
for further mathematical studies. 

In any society we actually deal not with one object of satisfac- 
tion but with a large number of them. There are also, as remarked 
above, many different kinds of work, each having for a given in- 
dividual its own optimal value y*. The production of an object of 
satisfaction in general requires different types of work. Therefore 
the number n,, of types of work is in general greater than the num- 
ber no of objects of satisfaction. We still obtain an inequality 
somewhat similar to (55). The conditions for such inequality to 
hold may, however, express a number of relations between the now 
very large number of variables xi bias 2s ne a Th) 5 Yk (k = 
1,2,...,n,) and especially the a’s. Instead of a single value for 
a, an individual now may be characterized by an n, x no matrix 
with elements a;,. An element a;, gives the contribution of an 
amount y;, of work of the kth type to the production of an amount 2; 
of object of satisfaction of the ith type. It is likely to be an in- 
teresting and difficult mathematical problem to determine how 
nearly for given satisfaction functions S and s and for a given 
distribution function N (21, Lenses Lan Vik» Yrree* :Yn,) @ State 
of complete satisfaction of everybody can be reached. 


This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of the University of Chicago. 
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By the same method as in the previous paper by the author (Bull. Math. 
Biophysics, 23, 55-68, 1961), the mobile electron system is studied in 
tryptophan which has donated an electron to another molecule. The ob- 
tained net charge distribution of tryptophan" is compared with those of 
the semiquinones of riboflavin (FMN , FMNH and FMNH, ). From this 
comparison, a Suggestion is obtained for the structure of each of two 
charge transfer complexes between riboflavin and tryptophan in neutral 
and acid solution. 


From the observations of I. Isenberg and A. Szent-Gyorgyi (1958, 
1960), it is known that riboflavin can accept an electron from in- 
dole derivatives. These investigators have found that the longest 
wave length absorption peak of a solution of riboflavin and tryptophan 
(or other indole derivatives) is at 490 my in neutral medium, and at 570 
my in acid medium. These wave lengths are almost independent 
of the molecule which donates the electron. This indicates that 
the wave lengths of these absorption peaks are a property of ribo- 
flavin. It is important to notice that the wave length of the absorp- 
tion peak is dependent on the H* concentration of the medium. This 
indicates that one or both of the N atoms number 1 and 5 of FMN, 
where the H*’s are attached to the FMN (Karreman, 1961, herein- 
after referred to as I), are involved in the stabilization of the ribo- 
flavin” radical by the tryptophan* radical. The interaction of tryp- 
tophan (and in general, other indole derivatives) and riboflavin 
does not seem to be explained simply by the donor strength of the 
former (A. Szent-Gyorgyi, 1959). Therefore, it is interesting to 
compare the net mobile electronic contributions to the electrical 
charges of the atoms participating in the conjugation of the semi- 
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quinones FMN , FMNH and FMNH; (cf. I), with those of the tryp- 
tophan* radical. The mobile electron system of tryptophan is the 
same as that of indole. Therefore, we will consider in the re- 
mainder of this paper indole instead of tryptophan. These mobile 
net electronic contributions to the electrical charges of the atoms 
will be called here, for short, the net charges of the atoms. In the 
papers by B. Pullman and A. Pullman, (1958, 1959) the mobile 
electronic charges themselves are called ‘‘electrical charges.’’ 

The charges of the atoms in the indole* radical which arises 
from indole by donating an electron, as e.g., to riboflavin, is rep- 
resented in Figure 1. 

From the comparison of the molecular diagram in Figure 1 with 
those in Figures 5, 8 and 10 of I, the following features can be noticed. 

(1) In FMN , there are great negative charges on N atoms num- 
ber 1 and 5. In this relation it is important to notice that whereas 
the O atoms number 11 and 12 have also great negative charges, 
their influence is probably to a great extent neutralized by the 
great positive charges of C atoms number 2 and 4, and N atom 
number 3. 

(2) In FMNH, there is a great negative charge on N atom number 
5, whereas the magnitude of the negative charge on C atom number 
9 is great fora C atom. The influence of the two negative charges 
of the O atoms number 11 and 12 is, in the case of FMNH, even 
less than in that of FMN’, because of the additional great positive 
charge on N atom number 1. 


H H 
H C 0.148 C 0.115 
4 3 
A atic ae -0.055 
Bs ie H 
sC 0.080 A 0.005 
4 ae 0.071 N 0.525 
H H 


FIGURE 1: Net contributions of the delocalized mobile electrons to the 
electrical charges of the atoms in tryptophan*. 
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(3) In FMNH3, there are great positive charges on N atoms num- 
ber 1 and 5, 

(4) In indole*, N atom number 1 has the greatest positive charge 
of all atoms, whereas C atom number 4 has the next greatest posi- 
tive charge. 

(5) The N atom number 1, the C atoms number 13 and 14, and N 
atom number 5 in FMN , as well as the N atom number 5, C atoms 
number 15, 16, and 9 in FMNH, have approximately the same rela- 
tive positions as N atom number 1 and C atoms number 8, 9 and 4 
in indole*. This is the case because the corresponding bond lengths 
and bond angles involved have the same values in a first approx- 
imation. Also in this respect, the O atoms 11 and 12 compare less 
favorably with N atoms number 1 and 5 in FMN , or with N atom 
number 5 and C atom number 9 in FMNH. 

On the basis of these observations, the following structures are 
suggested for two different charge transfer complexes between the 
iso-alloxazine ring of riboflavin and indole. 

I. The FMN” —- Indole* Complex (Complex 1). Here the FMN™ and 
indole* are in parallel planes so that N atom number 1 in indole* 
is close to N atom number 1 of FMN’, and C atom number 4 in in- 
dole is close to N atom number 5 in FMN’. The first and third 
atoms, as well as the second and fourth, attract each other strongly. 

Il. The FMNH - Indole* Complex (Complez II). In this case the 
FMNH and indole” are in parallel planes, so that N atom number 1 
of indole* is close to N atom number 5 of FMNH, and C atom num- 
ber 4 of indole* is close to C atom number 9 of FMNH. 

It may be important in this connection to point out that in FMN~ 
and FMNH the N atom number 5 has a large share of the electron 
which has been accepted (see e.g. Fig. 8, Pullman and Pullman, 
1958). As a consequence, the negative charge of N atom number 
5 is much greater in FMN and in FMNH than in FMN. This great 
negative charge of N atom number 5 leads to a strong attraction 
between this atom in FMN~ and C atom number 4 in indole” in the 
case of complex I. In the case of complex II, the great negative 
charge of N atom number 5 is, in particular, strongly attracted by 
the great positive charge of N atom number 1 in indole*. Asa 
result, there will be a great stabilization of both these charge 
transfer complexes. 

Due to the electrostatic influence of the charges in indole* upon 
the electrons of FMN~ and FMNH, the electron affinities of the 

atoms of these two will change. Consequently, the wave lengths 
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of the absorption peaks will change. To determine the new absorp- 
tion peaks, we have considered in complex I, as a first approxima- 
tion, only the changes in the electronegativities of the N atoms 
number 1 and 5 of FMN. These are the atoms closest to the atoms 
with the greatest charges in the indole’ radical. We have sub- 
stituted the charges of the atoms involved and their distances in 
Coulomb’s law. Then, using the potential energies obtained, we 
have calculated the changes in the electron affinities of the N atoms 
number 1 and 5of FMN . By solving the new eigenvalue problem, we 
have: obtained the following results. When the planes of FMN and in- 
dole’ in complex I are at the distance of 34, FMN has its longest 
wave length absorption peak at 480 mu. The distance of 3A is a 
reasonable estimate for the following reasons. For the transfer of 
an electron from indole to riboflavin, a necessary condition is that 
the electronic wave functions of these compounds overlap. This 
overlap, which decreases at increasing distance, may be expected 
to be still sufficient at a distance of 3A between the planes of the 
compounds. On the other hand, that distance cannot be much 
smaller than 3A for steric reasons. 

Similarly, in complex II we have considered as a first approximation 
only the changes in the electron affinities of N atom number 5 and 
C atom number 9 in FMNH. These atoms are again closest to the 
atoms with the greatest charges in the indole’ radical. The results 
show that when the planes of FMNH and indole’ in the second com- 
plex are also at the distance of 34, FMNH has its longest wave length 
absorption peak at 590 mu. In a neutral solution of riboflavin and 
indole, the concentration of FMN” will be greater than that of 
FMNH and that of FMNH?, This is due to the low concentration of 
the H* ion and to the stabilization of the FMN” as a consequence 
of the strong electrostatic attraction of FMN~ by indole’. How- 
ever, as shown in I, the more acid the solution the higher will be 
the concentration of FMNH and that of FMNH? 9» compared with that 
of FMN~. In the presence of indole* in acid solution, FMNH will 
be stabilized due to the electrostatic attraction of FMNH by in- 
dole’. However, FMNH} cannot form a complex with the indole* 
because of the strong repulaten between these two radicals. The 
great positive charges of the N atoms number 1 and 5 in FMNH} 
strongly repel the great positive charge of the N atom pigs! 1 ant 
the positive charge of C atom number 4 in the indole’ radical. 
Therefore, a complex similar to complex I above will not be sta- 
ble. This may also be expected for a complex similar to complex 
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II. Moreover, other kinds of complexes will not be formed because 
of the strong overall en between FMNH> and indole*. There- 
fore, in acid solution FMNH; will not be bin hired’ whereas FMNH 
will be stabilized due to if attraction by indole’. Consequently, 
the concentration of the FMNH radical may be expected to be the 
greatest in acid, even strong acid solution. From these consid- 
erations and the above mentioned results of the calculations, the 
absorption peak may be expected to be at 480 mu in neutral solu- 
tion, and at 590 my in acid solution of riboflavin with indole. These 
calculated wave lengths of the absorption peaks are, within the 
accuracy of the method, equal to the experimentally obtained peak 
at 490 mp in neutral and 570 my in acid solution. This suggests 
that the peak at 490 mp is due to the absorption of FMN sta- 
bilized by indole”, and the peak at 570 my due to the absorption of 
FMNH stabilized by indole’. 

Because of the limitations of the method, especially in its ap- 
plication to the spectra of radicals, one should of course bear in 
mind the approximate nature of the agreement between calculations 
and experimental data. The main contribution of this work, based 
on quantum mechanical calculations and simple geometrical ob- 
servations, is that it may throw some light on the relative posi- 
tions of the two molecules in the charge transfer complexes con- 
sidered. Such positional relations, combined with the results of 
quantum mechanically obtained electric charges, are probably ap- 
plicable to many other charge transfer complexes. 

Later publications (in preparation) will deal with similar rela- 
tions which have been obtained for the structure of charge transfer 
complexes of several compounds. Other relative positions of the 
components of charge transfer complexes will then also be studied 
in more detail. 
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The supplantation of one of two closely similar competing species by 
the other (Volterra-Lotka principle) is studied as an example, according 
to Volterra’s general theory of population dynamics, of the decay of a 
biological association of an odd number 2n+ 1 of species into one with 
an even number 2n, The three-species association (n = 1) is worked out 
in detail, with a demonstration of how a Volterra predator-prey cycle gets 
gradually deformed—when a slightly superior predator is introduced— 
into another new-predator-prey cycle, spelling a steady eclipse of the 
original predator. When n is made large, the close competitors being em- 
bedded in an association of many other species, a statistical treatment of 
the supplantation process can be given through the author’s statistical- 
mechanical theory of Volterra’s dynamics. The result is a probability 
law, changing systematically as time goes on, for the chance that the 
successful competitor’s population has any given amplitude; explicitly 
time-dependent measures of mean amplitude and mean frequency of oscil- 
lation of all populations can then be worked out. Throughout, the simpli- 
fying assumption is made that the competitors differ only as regards in- 
trinsic rate of self-growth. Two things are accomplished by viewing the 
competition as a decay of 2n+ 1 into 2n: the competitors are not ab- 
stracted from the rest of the biological world, and their population varia- 
tions are always oscillatory (with long-term rises and falls of amplitude); 
this is in contrast to the original Volterra-Lotka analysis in which purely 
static, and therefore ecologically unrealistic, population levels of but 
two species eventuate. 


1. Introduction. The Volterra-Lotka principle—stating that two 
closely similar species will not both indefinitely be able to occupy 
essentially the same ecological niche, but that the slightly more 
“*successful’’ of the two will completely supplant the other eventu- 
ally—has had a prominent place in ecological and evolutionary 
theory. A recent review by Hardin (1960; rebuttal by Cole, 1960) 
shows how extensively the idea, already clearly appreciated by 
Darwin, has percolated through a century of biological thinking. 
Yet the mathematical formulation of the principle has apparently 
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been practically limited to the original rudimentary discussions of 
Volterra (1926), Lotka (1932), and Haldane (1924). In these, but 
two species—of population sizes N;, Ng—are considered which 
compete according to a kind of extended Verhulst-Pearl scheme 
(dN, /dt)/N, = pt — wi (aN1 + bNe), ordinarily leading to extinction 
of one species and dominance of the other at a static population 
level. 

We should like to remark in the present note upon a simple way 
of viewing the extinction-domination struggles of similar species 
that does not, as heretofore, abstract the struggle from the rest of 
the biological universe but rather explicitly recognizes it and af- 
fords some surveillance of both the strugglers and their milieu. 
The basis for this is Volterra’s (1932; also D’Ancona, 1954) gen- 
eral theory of interacting species, in particular that aspect of it— 
seemingly dubious and artificial at first sight—foretelling of the 
tendency of decay of a biological association of an odd number of 
species into one with an even number. First, this type of transition 
will be examined in all detail for the case of three species (two 
strugglers interacting with a third species), where the gradual as- 
cendancy of one struggler over the other to give asymptotically a 
stable (not static) two-species regime is brought out by a specific 
example. Then we attempt a discussion of the strugglers when em- 
bedded in a large association; here a complete analysis is out of 
the question, and we appeal instead to statistical-mechanical con- 
siderations of a sort previously discussed (Kerner, 1957, 1959); but 
now the statistical mechanics is basically of the non-equilibrium 
variety. 

The picture that comes forth is simply that of Volterra-type os- 
cillations, or fluctuations, of all the coupled species: the success- 
ful struggler has a secular rise in his amplitude of oscillation about 
his slowly rising average population level, the two rises proceed- 
ing to fixed limits irrespective of how small his original numbers 
may have been (it is assumed, naturally, that the population is not 
so small that its fundamentally stochastic behavior dominates); the 
less favored struggler experiences a gradual damping of his oscil- 
lations to extinction, no matter how numerous his kind was initially; 
the rest of the species in association practically, though not quite 
completely, are unmoved by the whole contest, fluctuating in popu- 
lation in their own ways through its long duration nearly as if the 
strugglers were indistinguishable and constituted a single species. 
In the course of the adiabatically slow alterations of the strug- 
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glers’ population amplitudes, their frequencies of oscillation re- 
main substantially invariable. 

Our use of Volterra’s dynamics for odd numbers of species, in 
Summary, will perhaps show a way to look at such ‘‘dissipative’’ 
characteristics of ecological systems as the rise and fall of species 
within one comprehensive scheme wherein important ‘‘conservative’’ 
elements belonging to the ‘‘even’’ dynamics still play their illumi- 
nating central role; the ad hoc addition of specifically dissipative 
mechanisms, like Verhulst-Pearl terms, with their unsatisfactory 
implications of ultimately completely static population levels, is 
avoided in favor of at least a positive response to what may be the 
starkest fact of macro-ecology, the unceasing see-sawing of popu- 
lation numbers beyond the happenstances of season and migration. 

2. Three-Species Interaction. Volterra’s dynamics describes the 
population variations of species in association through character- 
istic growth coefficients for each and interaction coefficients ex- 
pressing binary interactions of basically predator-prey type. Con- 
sider three species, of whom the first two are similar and non- 
interacting; the Volterra equations read 


. ‘ 
No = 61Ni+ — 1M Ns, 
B: 
° 1. 
N, = €2No + — Oo N2Ns, 
Ba 
; 1 
Ns = €3N3 + —(%13N3N1 + X23 Na No), 
Bs 
where the reciprocal character of interactions signifies 0;; = — jj. 


The similarity of 1 and 2 is to mean in general that €;, Bi, 3; are 
close in numerical size to €2, B2, X32- If 1, 2 are predators and 3 
a common prey, £1, €2 <0, €3 > 0, while 431, M32 > 0. Or, if 1, 2 
are similar prey for the common predator 3, these inequalities are 
reversed. This is to say that we are placing 1 and 3 and also 2 
and 3 in the classical Volterra pattern, which would show, for either 
pair considered separately, the well-known periodic cycle of popu- 
lation oscillation visualizable as a closed loop in the N;, Vg plane. 

The near-equality of vital coefficients of 1 and 2 is unfortunately 
difficult to exploit in a comprehensive way; however, when only 
€, and €9 differ, corresponding still to interesting, if not alto- 
gether general, biological conditions, the analysis is very simple. 
We shall limit our further discussion to this case. 
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We have now at once, from the first pair of equations, 
Ni Na 
Ves 
and so N;=RNoe&* or Ne = pNye © (p = No (0)/Ni (0) = Ra2)s 
Thence we can study the N;, V3 system, 
Ny = €1N1 + &1Ni Na, 
Nz = €3Ng — %3 N1N3 (1+ pe ©), 


or the No, Ns system, 
No = €gNo + aoNoNs, 
Ns = &3Ng3 - A3NaNa(1+ Re) 


(here 1 = 2 = 31/B1 = 432/B2; X3 = 431/B3 = X32/B3 = — 413/B3). 
To keep a definite illustrative picture in mind, we can, for instance, 
think of the predacious Nog and its prey Ng alone occupying the 
field for ¢<0 and executing a Volterra cycle (€2 <0, &2 > 0; 
€3 > 0, &3 > 0); at = 0, species Ni, the brother-competitor to Ne, 
turns up, say a little better equipped (€; < 0 also, but €; > €2 by 
the small margin €> 0), though perhaps only in small numbers 
(p >> 1); then over the long term (¢ >> 1/€) exactly how does the 
original No-N3 cycle transform into the basically similar Ni-N3 
cycle that eventually reigns as Ne gets eclipsed? 

Let us introduce N; = Ni(1+ pe ©*)=N,U(t) into the Ny, N3 
system of equations, 


fi, (4 ~ 1) + 0Mims, 
Ns = &3Ng — X3NiNa, 

where €; = — €; > 0. Now, since 

es ae: 

1+ Re®’ 


we see that |u/u | always < €. By hypothesis, ¢€ is exceedingly 
small compared with €, or €3. The inequality is therefore strong, 


and we shall not miss very much at all if we simply drop U/U to 
give, in a close approximation, 


U 
= 


1 = — €1Ny + &1N1 No, 
E3 Ng ~ a3 Ni Nz. 


=a =. 


oo 
I 
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This states that Nas N3 go forward in the standard Volterra cycle, 


~ 


or aN se 
Met. 2 =. const: Ny 


goes 

As Ny oscillates strictly periodically, Ny oscillates with it but 
grows steadily from N, /(1 +p) to, Ni, while No = Ny P Cg INES 
falls, as it, too, fluctuates, from Ni(p/1 +p) to zero. The same 
description, of course, follows from the No, N3 system after de- 
fining No = No (1+ Re®*) and examining the Ny, Ng closed cycle. 

A slightly different method of analysis comes from noticing that, 
over any sufficiently small interval of time, 1+ pe~®* can be repre- 
sented as accurately as desired by a simple exponential Ae~**; and 
then, if Ny = N,Ae7At, the N;, N3 system becomes 


Nie TETAS ON Ne 5 


Gre. Nes NN G, 


Hl 


from which the usual Volterra integral follows for the duration of 
the small time interval. Thus for any of the intervals (0, ¢,), 
eae ie ip. ti p-.<, We can. place - 14 per@*kt*) 
1+p,e ©*(p, = pe ©**) and approximate it as (1+ p,)e “© with 
Lx < px/1+ px, in order that the latter approximation may make a 
little arc below 1+ p,e © that hits it at s = 0 and a little beyond 
at s =¢,4,; —t,- The ‘‘little’’ arc can, of course, be made quite 
long if € is small and p, large; for instance, the maximum error can 
be arranged to occur for s ~ €~' and then, if p, + 1-= 104, this 
error is some 0.1%, the arc extending from s = 0 tos ~ 6.5 ', and 
yu, being 0.99. The stepwise integration can be carried as far and 
as accurately as needed, but little is learned by so doing beyond 
what the approximation U/U = 0 already tells with enough clarity. 
In Figures 1-5 we illustrate what happens in a particular case. 
The Ni, Ns variables | are more conveniently replaced by the re- 
duced variables 7; = N1/Q1; ng =N3/Q3, where Q, = &3/H3 and 
Q3 = €1/%1 represent the ayerage levels about which M1; N3 Os- 
cillate; similarly, no = No/Qa with ae ake The variables in the 


figures are 


N ae kek No 
1 ig e = 01 1+Re 2 


which therefore measure the strugglers’ populations against the 
average | levels they would have if they weren’t in competition; 
these (0, and Q2) are also the levels at which each population 
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= 
n 3 


FIGURE 1. Variation of the populations of successful predator (n Ph. and 
prey (m3). The spiral continues to wind upward until closing in asymp- 
totically on another loop corresponding to ni’s, completely displacing the 
unsuccessful predator no (Fig. 3). The marked points show the closed 
71,ng3 cycle, and correspond to times indicated in Figure 2. The numeri- 
cal values €; = 1, €3 = 2, €= 0.05, p = 2.5 have been used. 


separately could be static, Ns being simultaneously static at 
Q3 = E1/o 1: 

3. Many-Species Interaction. Let us now immerse the strugglers 
N,, No in a large association and try to trace out something of the 
behavior of all species. We have 


Ni = EyNy + + 3 a51N5N1, 
By 


No = €2No + ioe Asa NsNo, 
Bo 
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Nn = ou Ne + FE chet NeMh ire 3a 4. St 1); 
and shall assume B; = B2, 51 = Osa, Xie = 0, but €; ~ &9, and 
also that Ni, NM and No, Ny separately are even associations pos- 
sessing stationary states. 

Again, No = pN,e “' in the same notation as before, and elimi- 
nating N2 to allow examination of the Ni, Nx, system, 


Ny = €,N; + gre te NaN, 


1 
: 1 “7 
Ne = && NE + By fou,(1+ pe ©*\N Ni; + &g,Ng Nz + eee 
k 
As previously, put Ni =NiU (t) to give 


(1) 


= th = 
Serene ea arian Na Ne ans Net .-- 3 
k 


O qr 27 
WJ 


FIGURE 2. Time scale for cycles shown in Figures 1, 3, 4, 5, @ is the 
angle turned through in Figures 1, 3 by the radius vector from the point 
(1,1) to a point on the closed cycles %1, m3, and 7%, n3, measured counter- 
clockwise from an axis parallelling the ng-axis. The time is in units 
reciprocal to those for the €’s. 
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Ns 


FIGURE 3. Variation of the populations of unsuccessful predator (n3) 
and prey (n3). The spiral continues to wind down until squashed flat, 
when species 1 has wholly supplanted 2. 


In a first approximation, U/U can again be discarded in compari- 
son with €,, the accuracy of this being the better, the smaller ¢ is 
and the longer the time elapsed. It will suffice to consider just 
€ > 0 and just N;, N;. We shall discuss briefly below the effect 
of retaining U/U. 

The reduced N;, Nx system is now an ordinary Volterra system. 
Oscillation of N;, Nz occurs about the stationary levels 71, qx de- 
fined by N; = 0 =, 

O= 6:81 + Ss gs, 

O= & B+ {Xun 7i + sega t---$, 
i.e., the levels 71, q, are just those—gqi, g,—of the pure N,, Ny 
association. 


The full range of our earlier statistical-mechanical theory can be 
brought into play at this point: we cannot tell in detail the motion 


of any single variable, but, through the use of Gibbs ensembles in 
the phase space, 
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| 2 


oo 
oe 
FIGURE 4. Variation of the populations of the competing predators. 


The dashed line indicates what their behavior would be if they were 
identical instead of slightly different. 


we can characterize the system as whole by parameters of thermo- 
dynamic type, and can deduce statistical properties of a curve of 
01, v; or of Ni, N; plotted as a function of time. 


The most useful ensemble is iba enonical one, 
mG 
6 
eo 3e. ’ 


where 


~ 


Pree a ee tle © — e%) = Gi + & Gy 
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FIGURE 5. Variation of all three species in the nj, ng, ng space. Pro- 
jections of the space curve on the coordinate planes give the curves of 
Figures 1, 3, 4. 


is a constant of motion of the reduced system, and T1 = 71f1 = 
q1 B1, Te =~ Gk Be. Here p gives the probability of finding the sys- 
tem in the configuration 31, 3; + 421; vx, Vk +@v_ When the sys- 
tem has run a long time—the variables v passing through many 
cycles of oscillation—and attained ‘‘thermodynamic’’ equilibrium; 
6 and ~ then play a role analogous to temperature and free energy 
for physical systems. As we have previously discussed, we can 
then use canonical averages as the means for reckoning the time- 
averages which interest us observationally. 
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Thus, for example, any average involving the variables Dige wd dass 
or N3, N4,..., of the biological environs of species 1 and 2 will 
make use only of the 2G, part of G; this part is the same as it 
would be in the pure N;, Nx system, i.e., the environs are oblivious 
to the struggle. This holds only to the extent of the validity of the 
above approximation |U/U| << €1; when U/U is not completely 
neglected, the canonical distribution and this conclusion are slightly 
modified. 

Now consider averages involving the struggler N;. A primary re- 
sult of our previous work, applicable here, can be written 


= N NiU 
e°} _ if = bai zat iL = NiU (t) = 
qi q1 
The bar signifies ensemble-average in the first instance, and then 
can be understood to be time-average, whence 


1=0. 


a, Hf 
N,U(t) = ;[ N,U(t)dt = 1), 


0 
for large T. Now U(t) for very small € can be reckoned a slowly 
varying function against N;(t); that is, when the characteristic de- 
cay time, 7p = €', is very long compared with any time of oscilla- 
tion 7; of N;(¢), it is meaningful to write 


‘| zie { YB ae 
= Nu@ae~ Us [ Meth 0 G\ Nes 
= 0 £ 0 


if it is understood that T; << T << To; therefore, 


q1 x 
Ree BN ett) 
in pee q; (2) 


We find in this way a ‘‘short’’ time average, or moving average, of 
N, that is sensibly time-dependent over very long times and use the 
special letter g* to indicate this special average. This state of af- 
fairs is shown schematically in Figure 6. 

It is clear now that we can write 0; = log (Ni/q}) and conduct 
canonical averages quite generally according to this same scheme, 
regarding them as moving averages over time intervals covering 
many cycles of motion but short on the scale of time over which 
q(t) varies appreciably. ‘‘Short’’ here may mean decades or cen- 
turies, ‘‘long’’ perhaps millennia. Altogether, g; can be counted as 
constant in the canonical average of any function of N;; and after 
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FIGURE 6. Schematic diagram illustrating the moving interval 7 in- 
volved in the canonical average, and the relatively long interval 79 of 
secular variation, together with the relatively short interval 7; of popu- 
lation oscillation. 


this average is performed, its secular variation found through its 
dependence on qj (¢). 

Let us notice here that substantially all our previous results on 
averages of amplitudes and frequencies of oscillation hold un- 
changed but for the replacement of g; by gj: For example, we have 


[Nac \ 5 Oe i eae : 
<> 1) = nen Or Ng gee 
q3 93) ey Bi (1+ pe a 

Ny N41 1 

—+ -1])log 4 = —, 

q1 1 = 


log “ =o (21)-log #1 (@ = digamma function) , 
i 

and the like. The characteristic parameter, 21, which is time- 
independent, again dominates the scene. What is happening is 
simply something effectively like a slow vertical stretching of the 
Ni-versus-time curve, so that the curve at an early time, when Ny 
is small and excursions from it are small too, has an appearance 
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through a cylindrical magnifying lens, laid along the time-axis, just 
like that seen at a later time, with the unaided eye, when Nix and 
the fluctuations around it have increased. 

The basic result is that the probability that 0, will be found in 
the interval 3,,%1+d%1 is 


qi Bi 
6 


P(®1)d01 ~ exp - (e°'~%1)d2, , 


or, in N-language, 


N N N Ni\ d(Ni/qt 
(or )a(=) ~ exp — q1 Bi (“ — log 3) eg) a) : 
6 q1 qi Ni/qi 


This is an invitation to write 2, = qi 81/0 as gi B1/6*, with 6* 
defined as 


pe 


0 

Puspes ‘ 

and to say (not quite accurately) that the ‘‘temperature’’ 6* effec- 
tive for N,’s fluctuations rises adiabatically from @/(1+ p) to the 
asymptotic limit 0 effective for the other species, as gj rises to its 
full height g; and N; becomes a full-fledged member of the as- 
sociation. But it is perhaps simplest of all to refer the values of 
N,; to the mean level g; prevailing in the pure N;, N;, association 
that is reached asymptotically, that is, calling n; the variable 
N1/q1, to write 


(1 + penratat 


er PT) 


—Et 
=e, (lip en sang Fix —s 
é . n1 dni, 


P(n1)dn, = 


to exhibit the secular drifting of N,’s probability distribution. In 

this view both 2, (0) and z, (@*), not 6* alone, are needed to de- 
scribe N;. The infusion of a deterministic element 1+ pe °’ into 
a probability law is the mark of the extra determinate knowledge 
that we started with, that No = pN,e~**. 

It is quite feasible to take any moving average involving N; and 
showing the slow U(t)-type time dependence, and explicitly to 
calculate further the ultra-long-time average. To do so would be to 

| wipe out valuable knowledge. What is much more interesting is to 
regard U(t) itself as having stochastic elements, namely, p and E, 
and to evaluate averages over these on top of the canonical moving 
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averages and to retain the explicit time dependence. For in any 
current view of the association we may not expect to know, for in- 
stance, the value of p =N.2(0)/N,(0) with any sharpness, since 
‘*Q’? here refers most likely to the dim past; this lack of knowledge 
had then best be recognized openly through some probability state- 
ment about p. This begins to pay a bit of extra respect to popula- 
tion numbers as the stochastic variables they really are and, in 
particular, is altogether compatible with the statistical-mechanical 
viewpoint that is contemplating all manner of initial population 
sizes as a matter of principle. 

We see, finally, that a way is open for viewing the complete as- 
sociation in a general state of secular flux that overlies the shorter- 
term fluxes. There is no basic limitation to just two close com- 
petitors in what we have done; many, or all, species can be cast 
into those roles of close competition as are engendered by muta- 
tions or other mechanisms, and the larger workings of the Volterra- 
Lotka principle, which is to say of certain aspects of biologic 
evolution of the whole association, brought under study. It is 
hoped to discuss this in a later report. 

4. Correction to the Canonical Distribution. To the extent that 
U/U has been neglected in comparison with €,, our discussion 
represents an approximation. It is of interest at least sketchily to 
try to proceed more exactly and to assess the effect of the U/U 
term. Granting still that this term is sufficiently slowly varying, 
we expect it to produce a slight non-stationarity of the density-in- 
phase. The methods for modifying Gibbs’s canonical distribution 
in such a case have been discussed previously (Kerner, 1959). 

Going back to equation (1), let us consider any current value of 
U/U and calculate currently stationary levels 


dis rh ere of A Wee ce 


By (e, - ae = EOsige 


B3& =- Oe3Is 


Write g* = q, + qs, where the qs’S are once again the static levels 
of the pure N;, Nx association satisfying By &% = -2Oskqs, so 
that 


Bi Git = ata, 
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{1} denoting the column vector (1,0,0,...,0); a the matrix of the 
:;; and {g*} the column vector (g1% eee Inverting a, iq} = 
B1(U/U) (a7!) {1}, or, because a7! is antisymmetric ((a~!),, = 0), 
91 = 0, 


, U 
% = Bi 7 Ae PALS ne 313) 4, 4..15 
Using the same phase space, 0, v;, as before, and defining 


C= hig '—Bi¢i% + 3 Bk Get” — Ba (Qk + Vk) Vk» 
k 


we can write the equations of motion as 


oe AY Ost an =: =e Ask OG’ cf 1, ee 
BiBs 9%, Psheo9Ve BiB. ov, 
Owing to g,’s explicit time dependence, G’ now is no constant of 
the motion. 

It is the basic premise of non-equilibrium thermodynamics that 
the slow variation of the parameter X (¢) = 8, U/U occurring in @’ 
results in a quasi-stationary Gibbs ensemble, 

Y-C’ 
Bee (1S NYS poll LB 


telling that the density-in-phase in near-equilibrium is near-canoni- 
cal (with the current value of @’ acting like a near-constant of 
motion). 

The presumed small deviation from currently canonical behavior 
is expressed in the correction RY, where B is, in general, 6-, X-, 
and v-dependent and is found from Liouville’s equation upon as- 
suming that derivatives and squares of X are negligible, giving 


Leficay! .ditt% wo caged 

cal gies ay Sk De a US| ~~ 4 

6 \dX ax OV; Ov1 
Since po is normalized, the po-average of B must vanish. Our pres- 
ent X(t) already involves the derivative of a slowly varying func- 


tion, so the neglect of X and X? is strongly justified. 
The characteristic equation for B reads, after computation of 


dW/0X, 0G’/0X, 
A.B, 0G’ OB 
pea said ee Pt DD Yer ’ 
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where 2/ = (8,9, + XA,)/0, @, = BrI/9, Yer = %sr/Bs B, and where 
the tilde over v; has been diopped: 

This can be dealt with simply in the low-6 limit, when Gel0-= 
wy, + (@,— &;)U~% + 5 z,v;, and po is effectively Gaussian, for 
then the B equation can be written 


iN 


A, OB 
ae re - XA,) -)); YsrBs (9s Us — As Xx) 00, 


r=3 


It suffices here to put 
B =%0b,B,(q,v,—- X Ar) + const. 
with 


1- bs 
a Yor tr bp = — oe in mil 1) 
(5,1 = Kronecker delta symbol), i.e., 


bie 


A, 
ee = (1 a 551). 
Ts 


Finally, arranging for B (po-average) = 0, 


e 6, BrQr(v, — Ur) 


- SOL re 20-841), 


where v, also means po-average. 

It is possible to advance to higher 6 by a similar procedure, first 
writing the exponential in dG’/dv, as a truncated power series, and 
expanding B as a multiple power series in the v’s. This unfortu- 
nately rapidly becomes unwieldy and opaque. We must postpone a 
fuller discussion of non-equilibrium until the analysis here can be 
improved. 


B 
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A mechanism is suggested in which a postulated substance can enter 
a Synaptic region only during the simultaneous action of a conditioned 
and an unconditioned stimulus. If this substance remains and produces 
a physico-chemical change, then a consequence is that the degree of 
learning will depend on the spacing of conditioning trials and the time 
constant of the process will be related to the rate of diffusion or trans- 
port of the substance along the axon. 


While neural circuits may very well play a role in conditioning 
and thus form a convenient basis for learning theories (Rashevsky, 
1961; Landahl, 1941, 1943), evidence requires that there also be a 
more permanent basis for this phenomenon (Libet and Gerard, 
1938). It is the purpose here to consider some consequences of a 
simple model for conditioning which involves some kind of struc- 
tural change (cf. Hebb, 1949; Milner, 1960). 

Consider a neuron A upon whose cell body terminate two sets of 
end-bulbs, one set arising from neurons which are acted upon by an 
unconditioned stimulus S,, the other arising from neurons stimu- 
lated by a stimulus S,. Let s, and 8s, respectively denote the 
synaptic connections of a representative ending from a neuron of 
the first group and from one of the second group. By definition, 
the stimulus S, results in the action of the neuron A. We shall 
suppose that the neuron A is of the type in which the cell body 
membrane discharges or ‘‘fires’’ except in those parts of its surface 
insulated by synaptic contacts which are not active (Landahl, 


*This research was supported by the United States Air Force through 
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1961). Prior to any conditioning, on the other hand, the condi- 
tioned stimulus S, cannot by itself produce an action potential in 
neuron A so that an impulse coming to neuron A from S- dies out 
before reaching the region sc. Only when S, and S, act together, 
so that impulses arrive simultaneously at A, does the membrane of 
A in the region at s, become excited. 

Suppose that there is a substance (or group of substances) in the 
end bulb at s. which is able to penetrate into the synaptic inter- 
face sc when the membrane is electrically active, but otherwise 
cannot. Such a penetration might be the result of permeability 
changes associated with the action spike or it might be the result 
of the membrane current flow. A large number of such firings 
could result in an accumulation of the substance in the synaptic 
interface. If an accumulation of such a substance can cause an 
increase in synaptic area or a change in the electrical properties 
of the synapse, then a sufficiently large number of repetitions of 
S, and S, could make it possible for the stimulus S- by itself to 
fire A, and this effect could be relatively permanent. If reverber- 
ating circuits which are relatively short-lived are involved, then S, 
and Su need no longer be stimulated at exactly the same time and 
the number of repetitions of S, and S, can be greatly reduced. 

Thus we suppose that the substance moves from the axon at an 
approximately constant rate p into the end-bulb, where it is con- 
sumed in proportion to the amount, S, which is there. For this pur- 
pose we neglect any possible volume changes of the end-bulb. In 
the steady state and in the absence of conditioning there will be 
an amount So = p/y, y being a proportionality constant. But when- 
ever S, and S, act simultaneously an additional loss of the sub- 
stance occurs, the amount lost being proportional to the amount of 
S present at the time ¢. If Q is proportional to the average fre- 
quency of simultaneous firings of S, and S, resulting from some 
conditioning situation, the average being taken over a relatively 
long interval of the order of magnitude of 1/y, then we may write 
the following equation for the amount S of the substance: 


dS 
WP - YS - QS. (1) 


The measure of the permanent change is the amount of S incor- 
porated into the membrane, i.e., the time integral of QS. The 
amount of ‘*permanent’’ learning will depend on the amount of the 
permanent change. We shall suppose that the change in threshold, 
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or the effective change in excitation, is proportional to this inte- 
gral. If we denote that change by L, then the measure of learning, 
L, will be given by 


t 
af QSdt. (2) 
0 


Note that L corresponds to A€, = bc, for example, in previous for- 
mulations (Landahl, 1941; 1943). 

In order to investigate the consequences of the model, consider 
the situation in which n conditioning trials are spaced equally over 
an interval of time 7. Unless T is very large the rate @ will be 
proportional to the rate 


nh=n/T, 


so that for each trial there is, on the average, the same number of 
simultaneous firings of S- and Sy. From (1), with Q = 67, 6 being 
a proportionality constant, we have 


oe Spe (7 tt E(t = e (7FQ)4), (3) 
y+@ 


Hence from (2) and (3) we obtain L for ¢ 2 T from the expression 
AO pn 


SEP wee 
See y(n6 + yT)? 


[(n6 + yT)yT +nO1-e"*%T)] (4) 
so that the extent of learning depends on the duration of the inter- 
val over which the trials are given as well as on the number of 
trials. 

It should be pointed out that in the above model, interest is 
focused on relatively long term effects. A somewhat similar result 
can be obtained for the effect of short time intervals if such short 
times are considered that the activity of the reverberating circuits 
has not ceased before the beginning of the next trial. 

From equation (4) it can be seen that for small enough values of 
n or for large enough values of 7 the value of L is the same and is 
proportional to the number n, i.e., 


L (0, T) = L(n,x) = 46 pn/y. (5) 
If we compare the extent of learning for a constant number of 


trials during varying intervals 7 with the amount given by (5) we 
obtain the relative efficiency given by the following expression: 
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L(x) One yp? (One yt 


ya (1 - Pah iS (6) 


Since the only parameters are 6 and y, which are scale factors 
for the variables n and 7 respectively, expression (6) can be illus- 
trated graphically. This has been done in Figure 1, in which the 
efficiency is plotted against y7, with 6n as a parameter. From the 
figure it can be seen that according to this model the effect of 
spacing trials is to increase efficiency and this increase is more 
pronounced if there are many trials. 

It is a necessary consequence of the above model that the output 
neuron must become activated in order that a permanent change is 
able to develop. While it is likely that the resulting stimulation of 
a ‘‘pleasantness center’’ would enhance the rate of conditioning, a 
situation which can be taken into account, we could have consid- 
ered the case in which activity of the last neuron, A, does not lead 
to an overt response. 

As an example consider two stimuli S, and S,’ which are unable 
to produce the response # normally produced by the unconditioned 
stimulus S,. Let the presentation of S- result in afferent impulses 
which move along a chain of neurons, represented by a single 
neuron Na, so that the activity reaches, but fails to excite, the ef- 
ferent neuron A. Similarly, let S$,’ result in impulses which arrive 
at neuron NV,’ which in turn sends impulses to the efferent neuron 
A without, however, being able to activate the latter. 


1.0 
gnz0 
> 
Se gnz! 
Ee 
u 6 
hT gnz9 
J 
2 
| oa 4 
s 
“st gnzlO 
ae 
(0) 
0 | 2 3 4 3 6 7 8 9 10 


DURATION VT 


FIGURE 1, 


LEARNING 163 


Let there be weak connections from S,’ to Na and vice versa, SO 
that S.’ can send impulses to A but cannot excite Na. Let the 
stimuli S, and S~ be presented simultaneously a number of times. 
If neuron Na and the synaptic connections from S,’ have the neces- 
Sary properties, then the neuron A will be excited if SY is there- 
after presented alone. Now let S, be conditioned to the response 
to S,. After this has occurred, then, without any previous simulta - 
neous presentations of S,’ and S,, the presentation of S¢ will re- 
sult in the unconditioned response. It should be pointed out that 
such transfer of conditioning could also occur as the result of the 
activity in neural circuits. However, to test for this possibility in 
a conditioning situation is difficult, since S, and S,’ must gener- 
ally be presented a number of times together before conditioning of 
S. to S, takes place and hence a considerable length of time may 
have to elapse before a test for transfer is made. If S, and S¢ are 
presented simultaneously after conditioning of S, and Su, then 
either the unconditioned stimulus is not presented so that extinc- 
tion will occur, or the response occurs at the same time as do S, 
and S’ with the result that direct conditioning will obscure the 
effect of transfer. 
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A RELATIONAL THEORY OF THE STRUCTURAL CHANGES 
INDUCED IN BIOLOGICAL SYSTEMS BY ALTERATIONS 
IN ENVIRONMENT* 


ROBERT ROSEN 
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It is shown that a wide variety of structural alterations in both the 
*‘metabolic’’ and ‘‘genetic’’ apparatus of (Jl, 2)-systems can result from 
specific changes in the environment of such systems. A number of spe- 
cific examples are investigated in order to demonstrate the scope of these 
alterations. Certain biological applications of this discussion are sug- 
gested, including a suggestion for a possible interpretation of the mitotic 
cycle. 


We have pointed out elsewhere (Rosen, 1959a) that if the map- 
pings occurring in the abstract block diagram of some (il, &)-system 
A are not constant, then it becomes possible for suitable altera- 
tions of the environment of the system to induce various changes in 
the structure of A. It is the purpose of the present note to investi- 
gate this possibility in greater detail, and point out some of the 
possible biological applications of such phenomena. We assume 
that the reader is familiar with our previous work in this direction 
(Rosen, 1958a, 1958b, 1959a) and we shall adopt the notation and 
terminology of that work without further comment. 

It is best to proceed by the study of particular examples. Ac- 
cordingly, let us begin by considering the simplest possible (Jl, ®)- 
system, which consists of a single mapping f: A — B and a map- 
ping ®;:B —. H(A,B). We assume that this system is placed in an 
environment consisting of an element aeA such that ®;(f(a)) = f; 
i.e., the system functions as an (Jl, ®)-system in a completely sta- 
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ble manner. We investigate the behavior of this system in a chang- 
ing environment by allowing a to vary in A. In order to simplify the 
discussion, we shall for the time being ignore any changes in the 
various lags of the system caused by the environmental change. 

We now suppose that the environment of the given system con- 
sists of an element a’cd, where a’ 4a. Let us dispose of the 
trivial possibilities first. To begin with, it can happen that 7 is 
not-defined on a’, in which case the environmental alteration has 
destroyed all metabolic activity of the system. We shall henceforth 
ignore this most drastic possibility, although we point out that this 
possibility exists in each of the more complex situations which we 
shall consider below. 

Next, we may have 


f(a’) = f(@), (1) 


or more generally, 


®;(f(a’)) =f. (2) 
In the situations specified by (1) and (2), itis seen that the altera- 
tion of environment (i.e., the replacement of a by a’) has resulted 
in no significant change in the metabolic activity of the system. 
However, we remark that from the point of view of fine structure, 
case (2) is distinguishable from (1) by the fact that, in situation 
(2), f(a’) is in general different from f(a). 

We obtain more interesting effects when it is the case that 


O; (f(a’)) =f" #f- 
There are now various further possibilities. We may have 
f(a’) =f(@), (3) 


or more generally, 


®; (f’(a’)) =f. (4) 
In these situations, there is no great change in the activity of the 
system; however, it will be seen that the structure of the system 
will undergo periodic changes. That is, for certain periods, the 
metabolic activity of the system will be carried out by the com- 
ponent represented by the mapping f, and at other times by the com- 
ponent represented by the mapping f’. 

Instead of (4), we may have 


®;(f’(a’)) =f’. (5) 
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In this situation, we see that the effect of the change of environ- 
ment from a to a’ has been to cause the original system with which 
we began to become changed into a different system, possessing a 
different “‘metabolic’’ apparatus but retaining the same ‘‘genetic’’ 
apparatus. Moreover, we observe that once such a structural altera- 
tion has occurred, a restoration of the original environment need 
not be capable of causing a reversion to the original form of the 
system. In fact, such a reversion will occur if f’(a) = f(a), or more 
generally if ®s(f’(a)) = f. As a biological manifestation of morpho- 
logical changes caused by a change in environment of the above 
type, we may cite the gross morphological alterations caused by 
changing the kind and amount of food material in cultures of vari- 
ous protozoa (cf. Kudo, 1954, p. 109, for illustrative examples and 
further references). 
Finally, it may happen that 


®;(f"(a’)) =f” (6) 
where f’’ is different from both f and f’. In this case, it will be 
seen that the effect of the environmental change is to cause the 
system to ‘‘hunt’’ through the set of mappings H (A,B); i.e., the 
system will synthesize a sequence of mappings /%, f’, f®,...f™,..., 
with the sequence terminating at a mapping f° such that either 


a. ; (f) (a’)) = fo 
or 
b. ; (f° (a’)) = fr, 1Sk<n. 


In situation (a), it will be seen that the behavior of the system ulti- 
mately becomes completely stable in form and activity. Therefore, 
(a) is a generalization of the behavior specified by (5) above. On 
the other hand, in situation (b), the changes in the structure and 
form of the system ultimately become periodic (in much the same 
way that the decimal expansion of a rational number ultimately be- 
comes periodic); the period for the cyclical alterations of structure 
is clearly no —k. If it should happen that no mapping fo satisfy- 
ing either (a) or (b) is produced in the course of the system’s ac- 
tivity, then the structure of the system is both unstable and aperi- 
odic. We observe that the situation specified by (b) is a generaliza- 
tion of the behavior exhibited by systems satisfying (4) above. 
Next, we must investigate the effects of the alteration of en- 
vironment on the process whereby the ‘‘genetic’’ material of the 
system (i.e., the mapping ®,) is replicated. In order to implement 
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this investigation, we shall utilize the schema for such replication 
which was described in earlier work (Rosen, 1959a). We have 
pointed out subsequently (Rosen, 1959b) that this schema is proba- 
bly not completely appropriate to represent all the complexities of 
the biological situation; nevertheless, we believe it to be suffi- 
ciently representative of such schemata to provide a useful illus- 
trative example in the present context. 

We recall that the schema in question involved the embedding of 
the various outputs arising from the metabolic activity of an (Il, %)- 
system into sets of mappings appropriate to construct replicas of 
the ‘‘genetic’’ apparatus of the system. In the particular system 
which we have been considering, this schema requires that we as- 
sociate with each element 6€B a mapping B,, where 


B,:H(A,B) — H(B,H(A,B)). 


It is readily verified that we always have B;(f) = ®s, where of 
course f and ®; are the mappings which comprise our (il, R)-system 
in its original environment, and 6 = f(a). Thus, under the initial 
environmental conditions, ®; is replicated perfectly. However, 
since our replication schema depends intimately both on the struc- 
ture of the corresponding ‘‘metabolic’’ apparatus and on the en- 
vironmental outputs of that apparatus, it is to be expected that 
environmental alterations will be capable of eliciting aberrations in 
this replication process. 

This dependence of replication on structure makes it important 
for us to know precisely when, relative to the ‘‘metabolic’’ activity 
of the system, the replication process takes place. This problem 
has not yet been investigated, although it is our intention to return 
to the matter in a subsequent paper. For the present purpose, we 
shall henceforth assume that replication in the altered environment 
occurs before any environmentally induced alterations of structure, 
of the types specified by (3)-(6) above, have taken place. The de- 
tailed discussion of situations in which this is not the case is 
complex but not difficult and is left to the reader; we remark only 
that no essentially new situations are obtained by considering 
these situations. 

Now let us suppose once again that the environment of our sys- 
tem has been changed from a@ to a’. The environmental output of 
our system then changes from f(a) to f(a’), so that the induced rep- 
lication map is thereby changed from Bria) to B a7). There are 
then two possibilities: 


RELATIONAL THEORY OF BIOLOGICAL SYSTEMS 169 
Bray (f) = Of (7) 


or 

Bea (f) # Of. (8) 
Case (7) is trivial; no change in the structure of the ““genetic’’ ma- 
terial is elicited by the alteration of the environment. On the other 


hand, case (8) opens up a further set of possibilities, depending on 
the behavior of the mapping ©’ = By (a’)(f). It may be that 


®*(f(a’)) =f (9) 
O’ (f(a) =f" #f. (10) 


In the situation (9), we have the result that our replication schema 
has, in the changed environment, produced a ‘‘daughter’’ system 
which differs from the ‘‘parent’’ system in its ‘‘genetic’’ structure, 
but not in its metabolic structure. The situation (10) results in the 
production of a “‘daughter’’ system which differs from the ‘‘parent’’ 
system in both genetic and metabolic structure. The further be- 
havior of the daughter system arising in this case depends upon the 
relation between ®’ and f’. It can, for example, happen that 


Byca*y(f") = ®”. (11) 


In this case, it is seen that the structure of the ‘‘genetic’’ appara- 
tus in successive daughter systems will alter periodically, or 
**hunt,’? just as did the metabolic structure of the system in the 
situations described by (5) or (6) above. Moreover, just as in these 
situations, the various alterations in the structure of daughter sys- 
tems may have a finite period, or continue indefinitely. Perhaps 
the most instructive special case of this type of behavior occurs 
when the ‘‘period’’ of structural alteration is equal to two, and 


Byer ca’y(f") = OF - (12) 

In this case, the genetic structure repeats in successive genera- 
tions. This type of behavior is strongly reminiscent of the alterna- 
tion of generations which occurs in plants and certain invertebrates. 
We have thus far considered only the simplest possible type of 
(tl, 8)-system. However, it should be reasonably clear that basi- 
cally the same type of effects will occur in more complicated 
(St, 8)-systems as a result of environmental alteration. However, 
the greater richness of structure in the more complex (Jt, ®)-system 
will make the corresponding effects less simple. Since in general 
a complex (Jl, R)-system will require a large number of different 


or 
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sets to serve as environmental inputs, it is evident that an unlimited 
number of different effects of the above type can arise. It may 
happen that only certain components of the ‘‘metabolic’’ apparatus 
of the system will undergo the alteration in form we have described 
above, while others will remain unaffected; the same is true of the 
‘‘oenetic’’ apparatus of the system and of the induced mappings 
which function in our replication schema. Moreover, in the more 
complex systems, time lag alterations may play a prime role; we 
state without proof that it is possible for the structure of an (M, R)- 
system to be completely changed by means of comparatively simple 
changes in time lags arising from an altered environment. It is 
even possible for a connected system to break down into two or 
more fragments, each of which functions as an (Jl, )-system. We 
shall not now provide detailed examples of these various possibili- 
ties; the reader may find it of interest to undertake the construction 
of suitable examples for himself. 

Let us now examine briefly some of the biological applications 
of the foregoing discussion. We should point out, to begin with, 
that the specific effects of environmental alteration upon the meta- 
bolic and genetic structure of biological systems is a well-known 
phenomenon. It has been shown (cf. Benzer & Freese, 1958; Freese, 
1959) that specific environmental changes are capable of eliciting 
uniquely determined changes in the structure of the genome of ap- 
propriate organisms. In such circumstances the new environment is 
termed ‘‘mutagenic;’’ it is in this manner that the various effects 
described above are to be considered. Therefore our discussion 
has little to do with the periodic heredity-environment controversies 
which recur in the biological literature. 

As a further example of the possible importance of the above 
considerations for biological problems, let us consider the special 
case in which an (lt, R)-system operates by means of a feed-back 
mechanism. For simplicity, we again restrict the discussion to the 
simplest system {f, ®,}, bearing in mind that much the same argu- 
ment will apply to more complex systems. The feed-back mechanism 
in this simple case means that the set B of environmental outputs of 
the system is a subset (or contains a subset) of the set A of en- 
vironmental inputs to the system. In such a situation, all that is 
required is to start the system operating; once the system has com- 
menced to ‘‘metabolize,’’ it proceeds in effect to supply its own 
environment, and hence will actually proceed to regulate its own 
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structural alterations in accordance with the above discussion. 
To illustrate the possible significance of such a feed-back regu- 
lated system, we point out that perhaps the most familiar series of 
cyclical morphological alterations occurring in biology is the mito- 
tic cycle. The question as to what initiates the mitotic cycle has 
never been completely answered, but the most successful hypothe- 
Sis is that each cell, in the course of its metabolism, produces a 
substance that “‘triggers’’ the initiation of the cycle. It has, how- 
ever, remained an open question just how the profound morphologi- 
cal changes characteristic of the cycle are brought about. Although 
it is no doubt premature to assert that the mechanisms described 
above can provide detailed insights into this basic problem at this 
time, it should be clear that the situation seems to have much in 
common with a feedback-regulated (il, ®)-system capable of under- 
going cyclical morphological changes in altering environments. 

It is also to be hoped that the flexibility of the (Jl, R)-system in 
undergoing ‘‘morphological’’ changes in altering environments may 
be applied to problems analogous to those occurring in differentia- 
tion and carcinogenesis. However, these problems are intimately 
connected with the theory of the interaction of (il, ®)-systems, and 
it is hoped that a suitable discussion of these matters will be 
forthcoming in the near future. 
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The addition of a drug to a system which contains an active transport 
mechanism may produce a large number of effects on both the active trans- 
port mechanism and the permeability of the system. A general model of 
such a system is presented, and various possible results of the drug’s 
addition are discussed. Methods of experimentally separating these ef- 
fects on the active transport mechanism and the permeabilities of the 
system are then listed and analyzed. 


Introduction. When a drug is added to a system which contains 
an active transport mechanism, it may affect the active transport 
mechanism, the permeability properties of the system, or both. This 
paper is concerned with some of the various experimental techniques 
that may be used to ascertain where the drug actually affects the 
system and also the difficulties involved in any of these techniques. 
To do this, a very general model which incorporates the known 
salient features of most of the biological systems will be discussed 
first. Mechanisms for active transport and permeability will be re- 
viewed, and some of the ways in which a drug may act on these 
mechanisms within the context of the general model will then be 
discussed. In this discussion it will be seen that the drug may 
have many different actions which produce the same experimentally 
observed result. Finally, some experimental tests that can be used 
to distinguish between the effects of a drug on the permeability 
and/or the active transport mechanism wil! then be discussed. 
Many of these tests have appeared in the literature, but it is hoped 
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that, by discussing all the tests via this general approach to the 
problem, the reader will be able to decide which tests are most 
useful for him and what their limitations are for his system. It 
should be emphasized that we are here discussing tests to deter- 
mine whether a drug has altered the permeability, the active trans- 
port, or both, of a system, and not tests to determine how or where 
it does so. 

Model. The model we shall use is shown in Figure 1. The physi- 
cal rationale of this model is as follows. The system we are con- 
sidering is made up of many submembranes which are grouped to- 
gether to form a single membrane with solutions on either side of 
it. We are explicitly limiting our consideration to systems which 
(1) possess a single active transport mechanism for the particle 
species, A, under consideration; (2) vary only in the transverse di- 
rection; and (8) are in the non-excited resting state, i.e., we are 
not concerned in this paper with the effect of drugs on the excit- 
able characteristics of a system. 

The active transport mechanism, which will be called a ‘‘pump,”’ 
is assumed to lie completely within one of the submembranes. If it 
does not occupy the entire submembrane, there will be a portion of 
the submembrane on either side of it which possesses certain perme- 
abilities, denoted by Pp, and P,_. Moreover, A may be able to tra- 
verse this submembrane independently of the pump by means of a 
shunt permeability, P,. Thus the membrane may be divided into 
three portions: the submembrane in which the pump lies, that por- 
tion of the membrane to the left of this submembrane, and that por- 
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FIGURE 1. Model of active transport system. See text. 
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tion to the right. The last two portions will have permeabilities 
which will be denoted as P; and P,, respectively. It should be 
emphasized that the latter two regions need not be homogeneous. 
Lastly, the solutions on the left and right sides of the membrane 
will be denoted by 7, and r,, respectively. Many biological sys- 
tems are special cases of this model. For example, the erythrocyte 
appears to have a Na” and K* pump with a finite P, for each, while 
all the other permeabilities are effectively infinite. Other examples 
will be discussed throughout the article. 

We are considering what the system looks like to a single parti- 
cle species, which we have denoted as A. Thus the fact that P, 
may be close to zero for A does not mean that this portion of the 
membrane will have a small permeability for any other particle 
species. It may actually have a large permeability. (For example, 
in frog skin [Koefoed-Johnsen and Ussing, 1958], it appears that 
the portion of the membrane which has a low permeability to K* is 
highly permeable to Na*.) Also, the active transport mechanism 
need not pump A exclusively but may also pump other particles, 
either in conjunction with or in competition with A. Examples of 
this will be discussed later in this section. 

The permeabilities P,, P,, Pp), and P,, cannot be ignored in com- 
parison to P, in all cases. If the total membrane were very thin 
(less than a few hundred angstroms), these permeabilities might be 
ignored, but this is not true if the membrane consists of many sub- 
membranes or if it is a thick membrane. For example, in nerve 
fiber, if 7, refers to the external solution, the k* P; cannot be ig- 
nored in the analysis of the data (Frankenhaeuser and Hodgkin, 
1956), while in the frog skin, if 7 refers to the corium side of the 
skin, then the K* P; is smaller than its P,. 

We shall now discuss what the terms ‘‘pump’’ and ‘‘permeability”’ 
actually imply. By ‘‘pump’’ is meant that portion of the system 
which can transport A against its electrochemical gradient. (This 
definition is by no means correct for all applications, as has been 
extensively discussed in the literature, but it is as good a defini- 
tion as any yet offered.) The pump may work in a large number of 
different ways (cf. Patlak, 1956). The question of whether or not 
there is a pump in the system is frequently quite ambiguous. For 
example, if there is a solvent flow that carries A along with it, it 
is a matter of choice as to whether one calls this solvent flow a 
pump or not. This will be discussed further later in this section. 
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The pump may transport A alone, or it may transport other parti- 
cles along with A. Consider the case where A is an ionized parti- 
cle, the membrane is highly impermeable to charged particles, and 
the pump is a carrier model. The carrier may be charged opposite 
to A, the same as A, or neutral. If oppositely charged, then it will 
carry enough A so that the complex is neutral. However, for the 
carrier to return to its initial position without carrying the A back 
with it, and thus yielding no net transport, it must carry along some 
other particle with the same charge as A, so that their complex is 


also neutral. This model has been proposed for the transport of 
Na* and Kt in the frog skin, erythrocyte, nerve fiber, and other 
systems. We will, for the sake of definiteness, call it a “‘double 
pump.’’? It should be emphasized that the other charged particle 
need not be unique but may be one of a number of particles. The 
choice of which particle is used will depend on both the concentra- 
tion of the particles and the affinity of these particles for the car- 
rier sites. For example, H* might compete with both Na* and K* in 
the above examples, but it may not play an important role because 
of its low concentration. If the carrier has the same charge as A, 
or is neutral, then it will transport A and oppositely charged parti- 
cles simultaneously, such that the complex is neutral. We will call 
this mechanism an ‘‘ion-pair pump.’”’ If the pump transports only A, 
we will call this a ‘‘single pump.”’ 

The pump need not carry A exclusively but will probably trans- 
port other similar particles as well. Thus, if the concentrations are 
large enough, these particles will offer serious competition to the 
transport of A. This fact has been used to determine various char- 
acteristics of specific active transport mechanisms. (Cf. Christen- 
sen, 1959, 1960, for a review of the application of this method for 
amino acid transport and J. F. Sutcliffe, 1956, for an application to 
plant cells.) 

We shall now turn to a discussion of permeability. By **perme- 
ability’? is meant that portion of a submembrane which allows A to 
cross the submembrane without being actively transported. The 
most common type is that in which the A diffuses across the sub- 
membrane under the ‘‘force’’ of its own electrochemical gradient. 
The A may also diffuse across in conjunction with another parti- 
cle, as in an ion pair. (Cf. Brodie and Hogben, 1957, for a review 
of membrane permeability to undissociated molecules of weak elec- 
trolytes.) However, there are many other mechanisms for allowing 
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A to traverse the membrane without being pumped. In all the pos- 
sible models of a pump, some provision must be made such that the 
system supplies energy to the pump to transport A against its elec- 
trochemical gradient. If no energy is required, then the mechanism 
will no longer be an active transport mechanism but will be a perme- 
ability mechanism. Thus any mechanism which may be used as 
an active transport mechanism can also be used as a permeability 
mechanism, provided that the energy-yielding step can be and is 
omitted. We shall denote these mechanisms by the same names as 
those used above, except that the term ‘‘non-energetic’’ will be in- 
serted before the term ‘‘pump.’’ Some examples of non-energetic 
pumps are facilitated-diffusion, as defined by Danielli (1954) and 
further expanded by Bowyer (1957), and exchange-diffusion (Ussing, 
1949a), which is a special case of a double non-energetic pump in 
which the carrier transports A instead of a second particle back to 
the first side. Also, if the pump ‘“‘slips,’? so that for a carrier 
model, for example, the carrier and carrier-complex can at times 
cross the membrane without passing through the energetic steps, 
this portion of the pump’s action would be equivalent to a perme- 
ability. 

In the double non-energetic pump, if the concentration of one of 
the transported particles is increased, it will not only increase the 
rate at which the other particle is transported but will also allow 
the other particle to be transported against its electrochemical 
gradient. (This mechanism has been quantitatively discussed by 
Rosenberg and Wilbrandt, 1957.) Similar arguments would also 
apply to the ion-pair non-energetic pump and to the diffusion of ion 
pairs. However, the decision as to whether this mechanism should 
be called a ‘‘permeability’’ for the other particle, or a ‘‘pump,’’ is 
arbitrary. This, then, is another example of the difficulty in mak- 
ing a decision as to whether a mechanism is a pump or not. In 
terms of the definition of active transport given above, A is ac- 
tively transported in both the cases discussed above, but the ques- 
tion as to whether this decision is a useful one is an important 
consideration. We shall leave it to the reader to make whatever de- 
cision he wishes. 

A mechanism for permeability may also have an energy input. For 
example, energy may be constantly expended by the system, in 
order to keep pores across a membrane opened. This will not trans- 
port A against its electrochemical gradient, but, in order for A to 
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diffuse across the membrane, energy must be provided by the 
system. 

Effect of Drugs. We shall now discuss some of the possible ef- 
fects that the addition of a drug may have on a system. This will 
be done not only to show the large number of possible results but 
also to emphasize the difficulty in attempting to ascertain what the 
effect is for any given situation. It should be remembered that a 
drug may produce any or all of these effects. 

The drug may affect the permeability in a number of ways. For 
the type of permeability which allows A to diffuse across the mem- 
brane, the drug may alter, among other things, (1) the structure of 
the membrane (for example, the size of the pore in membranes with 
pores); (2) the partition coefficient inside the membranes (so that 
the concentration of A inside the membrane will change); and (3) 
the fixed charge inside the membrane (which will alter the concen- 
tration of charged A’s). (In all cases of drug action, we shall not 
discuss any specific chemical mechanism of how the drug actually 
accomplishes its effect.) If the size of the pore is controlled by an 
energetic reaction, then the altering of this reaction will also af- 
fect the permeability. If the permeability is partially or totally via 
a non-energetic pump, then the drug may affect this pump in a way 
similar to the way it affects energetic pumps, as discussed below. 
The drug may also react with A in such a manner that the complex 
diffuses across the system. That is, the drug itself may act as a 
non-energetic carrier, and thus the permeability of the system to A 
would appear to increase. If the permeability of the system to A is 
very low, then even a low concentration of the drug could produce 
an appreciable change in the apparent permeability. Finally, the 
drug may react with fixed charges in the solutions on either side of 
the membrane or with any component that binds A. 

When a drug is added, it may not be able to permeate the entire 
membrane and therefore will affect only some parts of the mem- 
brane and not others. That is, its action on the permeabilities may 
be different if it is added on opposite sides of the membrane. The 
Same argument applies to the effect of the drug on the pump (Bowyer 
and Widdas, 1957). 

The drug may perform a large number of different actions on the 
pump. It may alter the rate at which the pump works by affecting 
the forward and/or backward rates of reaction of the pump. This 
may be done either while altering the equilibrium constant of the 
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pump, by affecting some of the reactants involved or while not al- 
tering the equilibrium constant, either by affecting only enzymes 
involved in the reaction or by altering the diffusion constant within 
a carrier model pump. If there is more than one energetic pathway 
through the pump, the drug could also act by affecting the different 
pathways and could thus actually cause a change in the velocity 
and equilibrium constant by decreasing theireaction rate in the usu- 
ally used pathway while increasing the rate in a different pathway. 

The drug may also convert the pump into a permeability. For 
example, if the drug can ‘‘uncouple’’ the energetic portion of the 
active transport mechanism, the resultant transport will be pas- 
sive. Practically, this will be applicable only to systems in which 
the energy steps are not involved in making normally slow forward 
reactions go faster, since, if these reactions are stopped, it will 
experimentally appear as if the pump had stopped completely. That 
is, the energetic reactions should be involved, for the carrier model, 
with ‘‘binding’’ or transforming the carrier so that A cannot be 
transported from the higher to the lower electrochemical side or, in 
a double pump, with putting the other particle species onto the car- 
rier to a greater extent than it would go on passively. 

A drug may have an indirect action on a system. For example, 
in a double pump, if the drug affects the rate at which one of the 
two substances is transported, it will most probably affect the rate 
at which the other substance is transported. This is not surprising 
because, from the viewpoint of A, the other particle species which 
is pumped in conjunction with A may actually be considered as 
part of the pump. Thus any action which affects the other sub- 
stance acts as if the pump for A had been altered. Another more 
extreme example is the situation in which a second particle spe- 
cies, although not pumped, plays a role in the reactions involved 
in A’s pump. The drug may affect this second particle species 
by altering the membrane permeability to it and therefore alter- 
ing its concentration, or, even more indirectly, the drug, by al- 
tering the permeability to a third charged species, may alter the 
electric potential and thus may alter the concentration of the sec- 
ond species if it is charged. This will alter the pump rate, and 
therefore it will appear as if the drug acted directly on the pump. 
A final example is where the drug alters the concentration of a 
second particle species, which will cause an osmotic imbalance, 
so that the resultant water flow transports A. 
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The drug may have many different effects on the unidirectional 
flow of tracers of A. Since it may alter the permeability and the 
pump independently, the drug can cause the influx and outflux inde- 
pendently to increase, decrease, or remain constant. However, it 
can also produce these results by its action on the pump alone. 
For example, if the net effect of the drug’s action is to decrease 
the net flow of A through the pump from the lower to the higher 
electrochemical region (or even reverse this flow), it can at the 
same time actually increase the tracer flow in this direction. An 
example of this, in a carrier mechanism, would be the case of a 
drug which depresses an energetic step involved in both freeing A 
from the carrier on the higher electrochemical side and at the same 
time ‘‘binding”’ the carrier site so that no A could attach itself to 
the carrier on this side of the membrane. If the free carrier diffused 
back at a much slower rate than the carrier-A complex, then, when 
the drug inhibited this reaction, tracer A on the higher electro- 
chemical side could now combine with the carrier, and, because of 
the larger diffusion rate of this complex over that of the free car- 
rier, there would be more carrier available to combine with the 
tracer A on the lower electrochemical side. Thus both tracer flows 
would increase. Other examples can be found in the model dis- 
cussed by C. Patlak (1957, eq. 8), if the drug is assumed to alter 
various reaction-rate constants. 

Since a drug can cause a large number of possible effects, it is 
extremely difficult to devise a single test which can discriminate 
between all these results. Of course, if the specific character- 
istics of a particular system are known, many experiments which 
rely on the peculiar properties of the system may be devised to as- 
certain the drug action. However, this approach will not be dis- 
cussed in this paper. 

A single measurement (for a single set of conditions) may not be 
sufficient to ascertain the drug action and may not even tell whether 
the drug had any action. For example, if the system is in pseudo- 
equilibrium (i.e., there is no net flow of any particle species across 
the membrane) and the drug alters P; and/or P,, there will be no ob- 
served change in the concentrations and potential across the 
membrane. 

A difficulty in deciding whether the drug even has an effect on 
the system arises from the presence of the rate-limiting steps in 
both the pump and the permeabilities. For example, if P; is large 
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compared to P, and Pp,» then if P; is increased further, no effect 
would be noted. If a drug should cause a decrease in the rate of a 
non-limiting step, its effect might not be seen at low concentra- 
tions of the drug, whereas at higher concentrations of the drug this 
step might then become rate-limiting and therefore observable. 

Another difficulty arises from the fact that there is experimental 
error in all experiments. Thus, if the drug makes P, very large (of 
the order of the permeability of an equivalent thickness of solu- 
tion), if P; and P, are much less than P,, and if the pump’s rate is 
very small, then it would be experimentally impossible to determine 
whether the drug affected the pump or not. This is because the 
flow through the pump, for any concentration level of A, would be 
less than the error in measuring the net flow through the perme- 
ability. Also, decreasing the amount of drug added, so that P, is 
not increased by such a large amount, may not help the situation, 
since the lower amount of drug could have had a smaller effect on 
the pump, and thus a null result, for the effect of the drug on the 
pump would not indicate anything. Another example is for the drug 
affecting an enzyme within the pump. This will not stop the opera- 
tion of the pump but will drastically decrease its rate. However, 
this action may appear experimentally as a complete stoppage of 
the pump. (That is, decreasing the rate of the pump by a factor of 
10 may appear experimentally as if the pump were completely 
eliminated.) 

Tests. We shall now turn our attention to some tests which may 
be applied to a given system to determine whether the drug affected 
the permeability and/or the pump. As stated before, we shall not 
discuss tests that are applicable only to very special situations 
but shall deal with those tests that may be useful for a number of 
different systems. The tests will be discussed in the following 
order: those having to do with (1) tracers; (2) time; (3) analogous 
responses; and (4) special attributes of the pump. We shall not 
discuss tests whose only purpose is to determine whether a pump 
is present or not (such as that derived by Ussing, 1949) but will 
assume that they are known, as there is a large body of literature 
on this point. Also, they can tell only whether the drug’s action is 
to stop the pump completely or not and cannot tell the effect of the 
drug on the permeability. 

1. Tracers. A rather general test involving tracers has been 
presented by O. Jardetsky (1957). The approach is to add the drug 
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and then note how the tracer flux of A changes from left to right 
and from right to left. Then, for all possible combinations of flux 
change, the article lists what the drug action is. However, there 
are certain rather limiting assumptions, besides those mentioned by 
the author, implicit in this approach. These are (1) the actively 
transported particle, A, is a neutral molecule; (2) the pump is, and 
remains after the addition of the drug, an irreversible pump; and 
(3) Pi, Pr, Pp,» and Pp, are all infinite. 

The first assumption that A is neutral is straightforward and can 
be tested easily. If A is a non-neutral molecule, then a change in 
the electric potential across the membrane could affect both the 
pump and the permeabilities in an unknown manner. Thus this ap- 
proach cannot be used for charged particles. 

The second assumption, of the pump’s irreversibility, although 
not rigorously correct in any case, may be experimentally valid for 
all practical purposes. That is, the equilibrium for the active trans- 
port reaction may be shifted far to the right. However, since a drug 
may affect the pump by either increasing or decreasing its forward 
pump rate, and at the same time independently increase or decrease 
its backward pump rate, all the results listed by Jardetsky may be 
explained on the basis of the effects on the pump alone. For ex- 
ample, the drug may affect an exchange-diffusion reaction (Tosteson 
and Johnson, 1957) which could actually be part of the pump itself 
(Heinz and Walsh, 1958) and would then, by this test, be taken as 
an effect on the permeability. Also, since this approach is de- 
signed for a single concentration of the A, it will not detect the 
effect of a drug whose action on the pump is only apparent at a 
different value of A than that used in the experiment. 

The third assumption, that P, is the only important permeability, 
is, as discussed before, not always valid. For those cases where 
it is false, this test may give erroneous results. Three examples 
of this are as follows (in all the examples to be discussed here, it 
will be assumed that the pump transports the neutral A only from 
left to right): (1) If one of the permeabilities is changed and if P; is 
not too large, the concentration distribution about the pump will 
change, and hence the rate at which the active transport mechanism 
pumps will be altered. This test would erroneously detect this as 
the drug acting on the pump. (2) The pump is oriented so that it 
picks up particles on its left and transports them to its right. The 
particles picked up, however, in a tracer study will be tracers 
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which originated from both the left and right sides (the relative 
numbers depending on the proportion of these tracers at the left 
side of the pump). If P; is not infinite in comparison with the other 
permeabilities, the number of tracers from the right side will not be 
negligible. Thus, if the drug causes the pump to pump faster (an 
analogous argument holds for a decrease in pump rate), it will pump 
more (in absolute value) of the tracers that originated on the right 
side. Thus the number of tracer particles that cross the membrane 
from the right side will be decreased. (A slightly different version 
of this effect has been proposed by Kirschner, 1955, 1959.) How- 
ever, this test would predict that the permeability was also affected 
(specifically, that it decreased). (3) If the concentration on the 
left was greater than that on the right and if the pump was ‘‘shorted’’ 
out (i.e., P, =), there would be a certain diffusion rate across 
P; and P, and a certain concentration in the region between these 
permeabilities. Now, in the actual system (P, # «), if the P; and P, 
are small enough so that the active transport mechanism pumps at a 
faster rate than the diffusion flow across the membrane with the 
pump shorted out, the concentration at the right of P; must be less 
than the concentration would be if P, were infinite, while the con- 
centration to the left of P, would be higher (since the concentration 
difference across a permeability must increase for the flow to in- 
crease). Since the concentrations to the right of P; and to the left 
of P, are equal when the pump is shorted out, the concentration to 
the right of P; in this case is less than the concentration to the 
left of P.. Thus there will be a flow across P, from right to left. 
Hence, when P, is increased, the backward flow across P, will in- 
crease and therefore the net flow across the membrane will de- 
crease. For the case of the tracer flow from left to right across the 
membrane, the concentration of the tracer on the right side is zero, 
and therefore the boundary conditions of the above discussion are 
fulfilled. However, the rate at which the tracer is pumped from left 
to right depends on the relative concentration of tracer from the 
right side. This factor complicates the argument a bit, but itcan be 
shown that the same result applies. That is, if the pump is fast 
enough, when P, is increased, the tracer flow from left to right ts 
decreased. This conclusion is contrary to the conclusion which 
would be drawn from this test. 

2. Time. Tests involving time as one of the variables are not 
usually applicable for a general system because of the unknown 
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time required for the drug to reach and act on the pump and perme- 
abilities. However, there is a test which has to do with the length 
of time needed for the effect of the drug’s action to become notice- 
able which can be used in a somewhat qualitative manner and thus 
may be applicable in a fairly large number of cases. The test, for 
electrolytes, is to see how fast the drug affects the electric po- 
tential across the membrane (Shanes, 1958). Some of the necessary 
assumptions involved in this application are that the drug reacts 
instantly and completely with the system and that the electric po- 
tential is determined solely by the relative permeabilities and ionic 
gradients across the membrane. If this is so, then a rapid change 
of the potential when the drug is added implies that the perme- 
ability of the membrane to the ions surrounding the membrane has 
been changed. On the other hand, a slow change can mean either a 
slow change in the permeability or a change of the pump’s action. 
The latter could be tested further by maintaining, if possible, con- 
stant concentrations on either side of the membrane and seeing 
whether the potential changed or not. If not, this implies no perme- 
ability change. However, if there were a change, this would imply 
either a permeability change or that the pump itself had been modi- 
fied such that ions could flow through it. The effect of the drug on 
the permeability to A could be further tested by changing the con- 
centration of A (or removing it entirely, if this is possible) and ob- 
serving how the potential differed with and without the drug. 

If there were a barrier to the diffusion of the drug or if it reacted 
slowly with the system, then a slow change in the potential would 
not imply anything. On the other hand, a fast change would imply 
not only a permeability change but also that the drug’s action was 
rapid. The assumption that the electric potential is determined 
solely by the relative permeabilities and ionic gradients across the 
membrane is a rather important one. This assumption means that 
the pump does not play any role in determining the potential except 
insofar as it contributes ions to form the ionic gradient. Thus both 
the double pump and the ion-pair pump would fit this description, 
whereas a single pump would not. With a single pump, any change 
in the rate of transport would instantly result in a change in the 
electric potential across the membrane. Thus, if a drug which is 
known to affect only the rate of active transport is added to a 
membrane, and there is no change or only a slow change in the po- 
tential across the membrane, then, unless the drug acts very slowly 
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on the pump (such as by a diffusion barrier to the drug reaching the 
the pump or by the drug not being able to destroy the high-energy 
phosphate reserves of the system, although it stops glycolysis and 
respiration, cf. Maizels et al., 1958) or there is an exact equal and 
opposite effect on the permeability at the same time, the pump can- 
not transport a net charge across itself. (This is one of the argu- 
ments in favor of a Nat — K* double pump in nerve fiber, Hodgkin 
and Keynes, 1955.) However, if the amount of ion transported by 
the pump is small and the permeabilities to the ions large, a change 
in the pump’s rate by a drug may be reflected in only a very small 
change in the potential, which may not be experimentally observed. 
(The size of the capacitance of the membrane is, contrary to the 
statement of Hodgkin and Keynes, loc. cit., not important in this 
analysis except insofar as it helps to determine the time scale 
needed to discriminate between a single or a non-single pump.) 

3. Analogous Responses. In the ‘‘analogous response’’ tests 
the effect of the drug is found by analogy to some other phenomena. 
One application of this test is in determining whether the drug af- 
fects the permeabilities for a different species of particles, as, 
for example, one which is not actively transported. If so, it is 
probable that the drug may also affect the permeability to the species 
of particle under consideration, although, of course, it is obvious 
that the drug may be specific for the different species permeability. 
A simple version of this method, for electrolytes, is to see whether 
the electrical conductivity of the membrane is modified by the ac- 
tion of the drug. However, there are a large number of problems 
connected with this approach, some of which are as follows. 

If the pump is of the type which yields a net transport charge 
(i.e., a single pump), part of the conductivity measured would be 
due to the action of the pump. This portion of the electric conduc- 
tivity may not be negligible if the rate of the active transport mech- 
anism is similar to the rate of movement of all the ions via their 
permeabilities. To avoid this difficulty, assuming that the pump 
has a maximum pumping rate, large currents could be passed through 
the membrane so that the effect of the pump would be negligible. 
(If the membrane is an excitable one, this modification cannot, of 
course, be used.) 

Another difficulty involved in the electric conductivity measure- 
ment approach is that the conductivity is a function not only of the 
permeability but also of the concentrations of the ions. Thus, if 
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the drug alters only the pump, this may change the concentrations 
so that the conductivity is changed, and therefore it may appear 
that the drug also affected the permeability. If the concentrations 
could be held fixed after the addition of the drug, there would be 
no problem, but in many cases (frog skin, erythrocyte, muscle, etc.) 
the internal concentration of the ion cannot be controlled or may 
not even be known. Also, since the length of time needed for the 
drug to act on the pump may not be known, a measurement of the 
conductivity coincident with the addition of the drug, or shortly 
after, may not prove useful. That the change of concentration may 
be important can be seen from the following example. 

Consider a system in which a membrane surrounds an inner re- 
gion. Assume a double pump with the P, permeable to the uni- 
valent cation, B, which is transported from the outer to the inner 
solution, and relatively impermeable to the other univalent cation, 
A, which is transported from the inner to the outer solutions. Fur- 
ther, assume no impermeable ions inside or outside. Thus 4 will 
be concentrated outside; there will be an electric potential dif- 
ference which is positive outside to inside; B and any other posi- 
tive ion will be concentrated inside, because of the potential dif- 
ference; and the negative ions will be concentrated outside, be- 
cause of the potential difference. If the external solution is kept 
constant and if the drug’s sole effect is to cause the pump to stop, 
then, upon poisoning the pump, the system will move toward an 
equilibrium in which there is no potential difference, and the in- 
side concentration will be equal to the fixed outside concentration. 
Hence, in this state, the concentration of the positive ions will de- 
crease inside, and the concentration of the negative ions will in- 
crease inside. (We may ignore the change in A, since we have as- 
sumed that the permeability to B is much larger than that to A and 
therefore A has a negligable contribution to the conductivity.) 
Thus the conductivity due to the positive ions will decrease upon 
the addition of the drug, while that due to the negative ions will 
increase. When this model is evaluated quantitatively (see Ap- 
pendix), it is seen that, for values of the potential that are biologi- 
cally plausible, there can be a net increase or decrease of the con- 
ductivities (depending on the relative values of the permeabilities 
of the positive and negative ions) by as much as an order of 
magnitude. 

Another version of the analogous response test is to observe the 
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effect of a drug whose action is known. Then add a second drug, 
either alone or in conjunction with the first drug, and observe the 
resultant effect. If the effect is different, then it may be assumed 
that the second drug acts on a different portion of the system than 
the first drug does. (This method, in a slightly more limited form, 
has been presented by Shanes, 1956.) We may note that tempera- 
ture changes, especially low temperature, may be regarded as a 
drug. That is, if the Qio of the pump is much higher than the Q1o 
of the permeabilities, low temperature may effectively ‘‘stop’’ the 
pump and alter the permeabilities only slightly. There are certain 
pitfalls to this method, four of which are discussed below. 

First, and perhaps the major one, is that the action of the first 
drug may not be exactly what it is thought to be. For example, in 
the article quoted above, anoxia and IAA are used to poison the 
pump. However, in a later work (Shanes, 1957) it was reported that 
anoxia also affects the permeability. Further, since the perme- 
ability may arise from non-energetic pumps, a drug which affects an 
enzyme system may also affect the non-energetic pump and there- 
fore affect the permeability. A second factor to note is that, even 
though the action of the second drug is different from that of the 
first drug, it may still be acting on the same portion of the system. 
For example, the two drugs may be acting on different reactions in 
the pump, with differing consequences, or they may even be acting 
on the same reaction but with differing degrees of effectiveness. 
A third factor is that, since drugs may have different actions at 
different concentrations, any difference in action may be due to 
this concentration effect. Lastly, the two drugs together may in- 
teract with each other or act on different portions of the same sys- 
tem such that their net result is different from the action of each 
one separately. Also if they have the same action, it may appear, 
upon addition of the second drug, as if the first drug’s concentra- 
tion was increased. Thus the effect of the two drugs alone should 
always be tested for. 

4. Special Attributes of the Pump. Certain tests may be used 
which depend on specific attributes common to most pumps. We 
shall consider two of these attributes: saturation and equilibrium. 

The saturation method uses the fact that if the pump’s action 
arises via a series of chemical reactions (as in a carrier mechanism 
of in a gate-type mechanism, as discussed by Patlak, 1957, and 
_ Mitchell, 1957), the rate of the reaction will be independent of the 
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concentration of A for sufficiently high concentrations of A. The 
test then is to increase the concentration of A on one side of the 
membrane until the pump is saturated, and therefore, when this con- 
centration is further increased, the increased flow will be due to 
the permeability of the membrane. Because of the unknown effect 
of the electric potential on the permeability and the pump for most 
systems, this test should be limited to neutral A’s, although it can 
be applied to systems in which the electric potential is small. 
Also, since non-energetic pumps can also be saturated, this test 
is good only for permeabilities which arise solely through diffusion 
pathways for the A, i.e., those permeabilities that obey Fick’s law. 
For these cases, the slope of the curve of flow versus concentra- 
tion for concentrations greater than saturation concentration will 
be the permeability of the membrane. Also, by use of this knowl- 
edge, the kinetic characteristics of the pump can be calculated 
from this curve. This method has been applied by E. Heinz and 
H. A. Mariani (1957) to the transport of glycine in Ehrlich mouse 
ascites carcinoma cells and by T. Shaw (1955) and I. Glynn (1956) 
to the transport of K* in erythrocytes. If a drug affects the perme- 
abilities P;, P,, and/or P, or affects the rate of pumping, it may be 
seen from a change in the shape of the curve. This approach was 
used by I. Glynn (1957) to test the effect of the cardiac glycoside, 
scillaren, on the passive permeabilities of the human red cell. 
However, the concentration of K* used was not large enough to 
yield unequivocal results. 

The equilibrium method uses the fact that if the pump’s action 
arises via a series of chemical reactions, then when these reac- 
tions are in equilibrium, the rate of transport through the pump is 
zero. To discuss this approach, we assume that there is a single 
transport pathway through the pump for the A and that the pump 
proceeds via a unique series of chemical reactions. Thus, letting 
My, Ne, and € represent constant coefficients, S, and P, represent 
substrates and products, respectively, involved in the reactions, 
and Aj, and A,, represent the concentration of A on the left and 
right side of the pump, respectively, then the over-all chemical re- 
action for the transport of A from the left to the right side of the 
pump is 


€ Ay, + miS1 + m2So + eee + mS = 


(1) 
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If represents the electrochemical potential of a substance, then 
the over-all free-energy change involved in transferring one mole of 
A across the pump is 
£F-y fo om 
av a) aa ar ie — i. g 
te ks if, cas "Se : 


r—O—~- 


Since the energy supplied by the system to transport one mole of 
the A from the left to the right is 


f Mi, zn 
W = a= Lees 
3 é Ys, ye ok (3) 
k=1 k=1 
therefore 
AF =p, — Hy) —-W. (4) 
Pp Pp 


If there is no shunt permeability, i.e., P, = 0, then when the chemi- 
cal reaction is at equilibrium, there will be no net flow of A across 
the pump and hence no net flow of A across the membranes and 
across the system. Thus the A in the left solution, A;,, will be 
in equilibrium with the A;,, and the A in the right solution, 4,,, 
will be in equilibrium with the A, . Therefore, 


(2.47) equilibrium a (#41, equilibrium? (5) 
(Ha 7 J ewdilibrivin - Ap egutiibrium ? 


Since AF = 0 at equilibrium, therefore equation (4) reduces to 
Wx (4, _ HA, equilibrium 


z 6 
wer in [eo HF] | - 
(A; ©) equilibrium 


where (A,.) and (A;,) are the activities of the A in the right and 
left solutions; z is the valence of the 4; V is the electric potential 
difference between the right and left sides; & is the gas constant; 
F is the Faraday constant; and T is the absolute temperature. (If 
the partition coefficient between the two solutions is not 1, the 
above equation may be modified in the obvious manner.) On the 
other hand, if P, does not equal zero, then when the system as a 
whole is in pseudo-equilibrium (i.e., there is no net transfer of any 
particle from one side to the other), the pump will not be in equi- 
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librium. That is, the A is traversing the P, from right to left, and 
the active transport mechanism is pumping A from left to right and 
thus at pseudo-equilibrium #4, 2 1 HA, Therefore, letting the sub- 
script ‘‘eq’’ denote the value of the variable in the pseudo-equilib- 
rium state, 


(41Jeq2 Aroq 224 (iA;.oq S (A,,)ea - 
Let 


sat a (A,*) 
Wap = (HA, in HAL dea ler ln (A, ) sy ; @) 
' eq 


Since A is being pumped from left to right, the reaction of equation 
(1) is proceeding from left to right, and therefore AF <0. Thus 


Woop = (H4,. = HA, Jeq $ (Ha, - HAL eq < We (8) 


6é 


The above result might be expected, since the pump must do 
tra’? work to overcome the leakage across it. 

An application of the above equations is as follows. If we fur- 
ther assume that the S’s and P’s are either electrically neutral or 
that those that have a charge do not move during the transport of 
the A’s, then the fs, and yp, will be independent of the electric 
potential. Assuming that the concentration of the S’s and P’s does 
not change during the course of an experiment where (A,.) and 
(A;,) and V are varied, we see from equation (3) that W will be con- 
stant. Thus, if P, = 0, equation (6) becomes 


Gas 
(A, a we constant. (9) 


eX- 


This relationship may then be tested by varying the A’s and the V, 
provided that the system is in pseudo-equilibrium. (In many sys- 
tems, the latter requirement may be experimentally difficult to 
satisfy.) If the relationship above is not constant, this implies 
that P, 4 0 (provided that the other assumptions are fulfilled). If 
the above relationship is constant, this does not, of course, prove 
that P,=0, but is consistent with this assumption. Hence, if 
equation (9) is satisfied, then the addition of a drug which then 
makes P, non-zero can be tested by this method. 

A further test is as follows. If equation (9) is satisfied (i.e., it 
is assumed that P, = 0), then if, after the addition of a drug, W 


app 
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increases, this can be explained only by the fact that the drug is 
affecting the equilibrium constant of one or more of the chemical 
reactions in the pump. The effect of the drug on P, could also be 
tested by seeing whether equation (9) was still valid. Note that if 
the drug merely affected the enzyme system involved in the active 
transport, but did not affect the equilibrium constants, this would 
not show up as an increase in W,,). 

For non-single pumps, some of the S’s and P’s represent the other 
particle species which are being transported by the pump. Thus 
the assumption that the charged S’s and P’s do not move would be 
violated, and equation (9) would have to be modified. For example, 
consider a double or ion-pair pump with the ion B, with valence 
2p, transported in conjunction with A. Then, with the same as- 
sumptions and with the P, of A and/or B equal to zero, we have, 
by an argument similar to that used to arrive at equation (9), 


A,. |78 [B,.|74 
= ts = constant. (10) 
Aj, eq B,, eq 


I wish to thank Professor Robert Macey for reading and discuss- 
ing this paper. 


APPENDIX 


We shall now quantify the discussion on the conductivity of the 
membrane with a double pump. To do this, we shall repeat our as- 
sumptions and add further ones as needed. Consider a single mem- 
brane, and, for simplicity, assume that P;, P,, Poy and Pp, are much 
larger than P, and that the ions in the solutions on either side of 
the membrane have a mobility which is much larger than their mo- 
bility in P,. Thus we may assume that the concentration and elec- 
tric potential are uniform in the entire region to the right and to the 
left of P,. (Therefore, the terms ‘‘P,’’ and ‘‘membrane’’ may be 
used interchangeably in this discussion.) Further, we will assume 
that P, arises solely through diffusion pathways for the ions and 
that the mobilities for the ions are constant within P,. We may 
group the ions into four classes: A, a univalent cation, which is 
pumped from left to right; B, a univalent cation which is pumped in 
conjunction with A via the double pump from right to left; other ca- 
tions which are not pumped; and anions which are not pumped. 
Then, if j, is the rate at which the &th ion diffuses across unit 
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area of the membrane; C; is the concentration of the &th ion; yo, is 
the standard chemical potential of the *th ion; y, is the activity 
coefficient of the &th ion; uz is the mobility in the membrane of the 
kth ion; V is the electric potential relative to the left solution; and 
2, is the valence of the &th ion, we have 


, dC. d ln yx d 10, dV 
=—_—? Sees Ss . 11 
ik mt rs + C; ler oe + Fa, a (11) 


Because of the difficulty in solving the above equations, we shall 
use the constant field assumption of Goldman (1943). This as- 
sumption is chosen not only because of its mathematically simpli- 
fying character but also because it has been found to be a good ap- 
proximation for thin membranes (Goldman, 1959). Letting AV repre- 
sent the electric potential across the membrane and 6 the thickness 
of the membrane (i.e., the thickness of P,) and assuming that both 
yx and yo, are constant, equation (11) becomes 


: aC; AV 
iE = — wm RT — + Oy Fay —1- (12) 
dz 5 
Assuming steady-state conditions, equation (12) integrates to 


aR igAV Ue Wa eae 


~~) 
5(1 —e RFA V/RT) : 


(18) 


In pseudo-equilibrium, if the ion is not pumped, j;, = 0. How- 
ever, if the ion is pumped, j, will be equal and opposite to the rate 
at which it is pumped. Since we are assuming a double pump for 
the univalent A and B, therefore j4 = —jg. Thus, by taking the 
sum of the A and B flow, we will get a total equal to zero. This 
fact can be used to find an interesting result. Adding the j4 and 
jp together from equation (11) without using the simplifying as- 
sumptions of the constant field or that the y’s and the zy’s are con- 
stant, integrating, and assuming only that y4 = yg, po 4 = Hog, and 
that both the y’s and the yo’s are the same in the two solutions on 
either side of the membrane, we have, in pseudo-equilibrium, 


RT w4Ca, + up Cp, 


AV = ph (14) 


u4Ca, + UB Ce, 


This equation has been derived previously for the double pump but 
only under the assumptions of a constant field and that the y’s and 
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Ho’S were constant (Hodgkin, 1957). However, we see that it is 
true under the much less restrictive assumptions listed above. 

We shall not enter into the argument of whether the conductivity 
is the ratio of the changes in the steady-state currents divided by 
the change in the electric potential, or the transient change in the 
current divided by the change in the electric potential. Since the 
latter is difficult to calculate and the former is easy to measure 
experimentally, we shall choose the former method. Thus the con- 
ductivity per unit area may be found by taking the negative of the 
derivative of the total current with respect to AV. Since the total 
current is simply the sum of 2, j,, then —2,0j,/dAV yields the 
contribution to the total conductivity from the sth ion. As the sys- 
tem is a non-linear one, the conductivities will not be constant as 
the AV is changed. Therefore, we shall choose to take the conduc- 
tivity measurement at the pseudo-equilibrium value. As stated be- 
fore, all of the j,’s will then be zero except for A and B, and, for 
this case, j4 + jg will be zero. Because of the mathematical simplifi- 
cation introduced by having the j’s equal to zero, we shall not 
calculate the conductivities of the A and B separately but will 
only consider them together. Therefore, from equation (18) and us- 
ing the fact that we are in pseudo-equilibrium, we have, when G;, is 
the partial conductivity of the éth ion, 


0 ix 2apF7uy, AV 
Ge = -—2k = = Cy, kK#A,B, 
GAV Re: (zp FAV/RT) r 
SRT(1-e ) (15) 
O (ja + iB) FAV 
eae aay is, SETA =CRAVIRTY, “A Ca, + us Cp,)+ 


If we add up the values for G, from equation (15) and make use of 
the fact that we are in pseudo-equilibrium, we will get the same re- 
sult as that found by Harris (1956, pp. 72-73). 

To apply the above results, consider a system, such as an eryth- 
rocyte, in which there is a membrane which surrounds an internal 
solution and is in turn surrounded by the external solution. When 
the pump is poisoned by a drug, the system will eventually reach 
the value given by a completely passive system. If we assume 
that there are no fixed or non-permeable charged particles in the in- 
ternal or external solution, then the resulting AV will be zero. For 
this case, letting C,, represent the concentration in the external 
solution, which does “not change during the course of the experi- 
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ment, the conductivities reduce to 


ae Fu 
G.(AVe=10) = ee k £ A,B, 
(16) 
F 
Gc asBy(AV = 0) = 5 (uw, Ca,+ UB Cp,)- 
Thus the ratio of the conductivities is 
Gy (AV) Zk FAV (17) 


G, (AV = 0) RTE (A — ee FAVRDY’ 


where the subscript & also includes (A + B). Choosing the biologi- 
cally plausible value for AV of 50 mV, the ratio of conductivities 
for univalent positive ions will be 2.27 and for univalent negative 
ions will be 0.33. If AV had been as high as 100 mV, the ratios 
would be 3.96 and 0.082, respectively. The assumption that ug > 
ua, although not made explicitly in the above derivation, is neces- 
sary so that AV will be positive. 


LITERATURE 


Bowyer, F. 1957. ‘‘The Kinetics of the Penetration of Nonelectrolytes 
into the Mammalian Erythrocyte.’’ /nt, Rev. Cytol., 6, 469-511. 
and W. F. Widdas. 1957. ‘*‘The Importance of the Reaction at 
Both Interfaces in Facilitated Transfer.’’ Jour. PAysiol., 136, 13-14P. 
Brodie, B. B. and C. A. M. Hogben. 1957. ‘*Some Physico-chemical Fac- 
tors in Drug Action.’’ Jour. Pharm. Pharmacol., 9, 345-380. 
Christensen, H. N. 1959. ‘‘Active Transport, with Special Reference to 
the Amino Acids.’’ Persp, in Biol. Med,, 2, 228-242. 
- 1960. ‘*‘Reactive Sites and Biological Transport Especially of 
Amino Acids.’’ Adv. in Prot. Chem., 15, 239-314. 
Danielli, J. F. 1954. ‘*‘The Present Position in the Field of Facilitated 
Diffusion and Selective Active Transport.’’ Colston Papers, 7, 1-14. 
Frankenhaeuser, B. and A. L. Hodgkin. 1956. ‘*The After-Effects of Im- 
pulses in the Giant Nerve Fibres of Loligo.’’ Jour, Physiol., 131, 
341-376. 

Glynn, I. M. 1956. ‘‘Sodium and Potassium Movements in Human Red 
Cells.’’ /bid., 134, 278-310. 

1957. ‘*The Action of Cardiac Glycosides on Sodium and Po- 
tassium Movements in Human Red Cells.’’ /bid,, 136, 148-173. 

Goldman, D. E. 19438. ‘‘Potential, Impedance, and Rectification in 
Membranes.’’ Jour. Gen, Physiol., 27, 37-60. 

1959. ‘‘Analog Computer Solution of Ion Flux Equations.’’ Ab- 
stracts of the Biophysical Society. 

Harris, E. J. 1956. Transport and Accumulation in Biological Systems. 
New York: Academic Press, Inc. 

Heinz, E. and H. A. Mariani. 1957. ‘‘Concentration Work and Energy 
Dissipation in Active Transport of Glycine into Carcinoma Cells.’? 
Jour, Biol. Chem., 228, 97-111. 


ACTIVE TRANSPORT 195 


and P. M. Walsh. 1958. ‘*Exchange Diffusion, Transport, and 
Intracellular Level of Amino Acids in Ehrlich Carcinoma Cells.’? Jour. 
Biol. Chem., 233, 1488-1493. 

Hodgkin, A. L. 1957. ‘‘Ionic Movements and Electrical Activity in Giant 
Nerve Fibres.’’ Proc. Roy. Soc., B., 148, 1-37. 

and R. D. Keynes. 1955. ‘‘Active Transport of Cations in Giant 
Axons from Sepia and Loligo.’’ Jour. Physiol., 128, 28-60. 

Jardetsky, O. 1957. ‘‘On the Distinction between the Effects of Agents 
on Active and Passive Transport of Ions.’’ Science, 125, 931-932. 

Kirschner, L. B. 1955. ‘‘On the Mechanism of Active Sodium Transport 
across the Frog Skin.’’ Jour. Cell. Comp. Physiol., 45, 61-88. 

1959. ‘‘The Interaction between Sodium Outflux and the Sodium 
Transport System in the Frog Skin.’’ /did., 53, 85-92. 

Koefoed-Johnsen, V., and H. H. Ussing. 1958. ‘‘The Nature of the Frog 
Skin Potential.’’ Acta, Physiol. Scand., 42, 298-308. 

Maizels, M., M. Remington, and R. Truscoe. 1958. ‘‘Metabolism and 
Sodium Transfer of Mouse Ascites Tumor Cells.’’ Jour. Physiol., 140, 
80-93. 

Mitchell, P. 1957. ‘‘A General Theory of Membrane Transport from 
Studies of Bacteria.’’ Nature, 180, 134-136. 

Patlak, C. S. 1956. ‘‘Contributions to the Theory of Active Transport.’’ 
Bull. Math. Biophysics, 18, 271-315. 

- 1957. ‘‘Contributions to the Theory of Active Transport: II 
The Gate Type Non-Carrier Mechanism and Generalizations Concern- 
ing Tracer Flow, Efficiency, and Measurement of Energy Expenditure.”’ 
Ibid,, 19, 209-235. 

Rosenberg, T. and W. Wilbrandt. 1957. ‘‘Uphill Transport Induced by 
Counterflow.’’ Jour. Gen. Physvol., 41, 289-296. 

Shanes, A. M. 1956. ‘‘Distinction between Effects on Metabolic Trans- 
port and Passive Transfer of Ions.’’ Science, 124, 724-725. 

———. 1957. ‘‘Ionic Transfer in a Vertebrate Nerve.’’ Metabolic As- 
pects of Transport across Cell Membranes, ed. Q. R. Murphy, pp. 127- 
150. Madison: Univ. Wisconsin Press. 

. 1958. ‘*Electrochemical Aspects of Physiological and Pharma- 
cological Action in Excitable Cells. Part I. The Resting Cell and Its 
Alteration by Extrinsic Factors.’’ Pharm. Fev., 10, 59-164. 

Shaw, T. I. 1955. ‘‘Potassium Movements in Washed Erythrocytes.” 
Jour, Physiol., 129, 464-475. 

Sutcliffe, J. F. 1956. ‘‘The Selective Absorption of Alkali Cations by 
Storage Tissues and Intact Barley Cells.’’ Potassium-Symposium 1956, 
pp. 1-11. Berne: International Potash Institute. 

Tosteson, D. C. and J. Johnson. 1957. ‘*The Coupling of Potassium 
Transport with Metabolism in Duck Red Cells: I. The Effect of Sodium 
Fluoride and Other Metabolic Inhibitors.’? Jour. Cell, Comp. Physiol., 
50, 169-183. “ 

Ussing, H. H. 1949a. ‘‘Transport of Ions across Cellular Membranes. 
Physiol. Rev., 29, 127-155. 

. 1949b. ‘The Distinction by Means of Tracers between Active 

Transport and Diffusion. The Transfer of Iodide across the Isolated 

Frog Skin.’’ Acta Physiol. Scand., 19, 43-56. 


RECEIVED 7-12-60 


SS SC 


ULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 28, 1961 


UBER DEN ERTRAGSVERLAUF BEI STEIGENDER 
NAHRSTOFFKONZENTRATION 
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In Nahrstoff-Steigerungsversuchen mit héheren Pflanzen oder mit Mikro- 
organismen wird in der Regel die Nahrstoffgabe 2 variiert, wobei diese 
meist vor oder bei Versuchsbeginn verabreicht wird. Man kann jedoch 
auch so vorgehen, dass man die Nahrstoffkonzentration u variiert, wobei 
die Konzentration u bei jeder einzelnen Variante wahrend der ganzen 
Wachstumszeit konstant gehalten wird. In der vorliegenden Arbeit wird 
untersucht, wie sich der Ertrag an Organismenmasse, y, bei steigender 
Nahrstoffkonzentration u verhalt. 

Bei Bakterien nimmt der ‘‘Ertragszuwachs’’ dy/du in jedem abge- 
schlossenen Konzentrationsbereich, in dem der Ertrag durch Erhdhung der 
Nahrstoffkonzentration u gefordert wird, entweder von Anfang an oder von 
einem wohlbestimmten spateren Zeitpunkt an mit wachsendem wu durchweg 
zu, sofern die Tochterindividuen immer wieder die gleichen Bedingungen 
vorfinden wie die Ausgangsindividuen. Ahnlich durften die Verhadltnisse 
bei anderen Mikroorganismen liegen. Auch bei hdoheren Pflanzen nimmt 
der ‘‘Ertragszuwachs’’ dy/du im aufsteigenden Ast dei Ertragskurve 
unter gewissen Voraussetzungen mit wachsendem u zu. Diese Ergebnisse 
erhalt man, wenn man die Abhangigkeit des Ertrages von der Zeit einerseits 
und von der Nahrstoffkonzentration andererseits im Zusammenhang 
betrachtet. 

Die gewonnenen Resultate widerlegen die weitverbreitete Ansicht, der 
‘‘Ertragszuwachs’’ nehme mit steigender Nahrstoffversorgung stets ab. 
Mithin ist das sogenannte ‘‘Gesetz vom abnehmenden Ertragszuwachs’’ 
kein allgemeingiltiges Naturgesetz. 


Einleitung. Die Frage, wie der Ertrag an Organismenmasse von 
der Nahrstoffversorgung abhangt, ist schon Gegenstand zahlreicher 
Untersuchungen gewesen. Trotzdem erscheint es notwendig, weitere 
Untersuchungen auf diesem Gebiet anzustellen. Im folgenden soll 
ein theoretischer Beitrag zu dieser Frage geleistet werden. Ehe 
mit den eigentlichen Untersuchungen begonnen wird, sei jedoch 
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kurz klargestellt, was unter N&hrstoffversorgung eigentlich zu 
verstehen ist. 

Bei der Durchfiihrung von N&hrstoff-Steigerungsversuchen mit 
hdheren Pflanzen oder mit Mikroorganismen wird gew6hnlich die 
Nahrstoffgabe 2 variiert, wobei diese in der Regel vor der Aussaat 
bzw. vor oder bei Versuchsbeginn verabreicht wird. In Stickstoff- 
Steigerungsversuchen beispielsweise stellt 2 die vor oder bei 
Versuchsbeginn eingebrachte Stickstoffmenge dar. Die N&hrstoff- 
gabe x wird etwa in Gramm oder in Milligramm je Gefass gemessen. 

Man kann Nahrstoff-Steigerungsversuche jedoch auch so durch- 
fiihren, dass man nicht die N&ahrstoffgabe z, sondern die N&hrstoff- 
konzentration wu variiert, wobei die Konzentration u bei jeder ein- 
zelnen Variante wahrend der ganzen Wachstumszeit konstant gehalten 
wird. Die Konzentration u kann man etwa in Milligramm oder in 
Millimol pro Liter Nahrldsung messen. 

Zwischen der N&hrstoffgabe x und der N&hrstoffkonzentration wu 
ist streng zu unterscheiden (Kaltofen, 1957a). Es ist nicht ohne 
weiteres angangig, x und uw einander gleich oder proportional zu 
setzen. Wird 2 variiert, so erhalt im Versuch jede Variante eine 
bestimmte N&hrstoffmenge, und zwar in der Regel vor oder bei 
Versuchsbeginn; hierbei ist die N&ahrstoffkonzentration in der 
Umgebung der resorbierenden Organismenoberflache im allgemeinen 
nicht standig die gleiche, da ja die Organismen N&ahrstoffe aufneh- 
men. Wird andererseits wu variiert, so wird durch geeignete Mass- 
nahmen dafiir Sorge getragen, dass sich die Nahrstoffkonzentration 
in der Umgebung der resorbierenden Organismenoberflache mit 
Fortschreiten der Zeit nicht andert; hierbei bleibt die Héhe der 
zugefiihrten Nahrstoffmenge ausser Betracht. 

Wenn man in den bisher durchgefihrten Na&hrstoff-Steigerungs- 
versuchen meist die Nahrstoffgabe z und nur selten die N&hrstoff- 
konzentration wu variiert hat, so hingt dies wohl nicht zuletzt damit 
zusammen, dass es technisch einfacher ist, vor oder bei Versuchs- 
beginn eine bestimmte Nahrstoffmenge einzubringen, als die Kon- 
zentration der Nahrlésung in der Umgebung der resorbierenden 
Organismenoberflache wahrend der ganzen Versuchsdauer konstant 
zu halten. Zweifellos gibt es aber Méglichkeiten zur Konstanthaltung 
der Konzentration. Was hdhere Pflanzen betrifft, so haben unter 
anderem W.-U. Behrens (1928) und K. Pirschle (1931) gezeigt, wie 
man die Ni&hrstoffkonzentration in der Umgebung der Wurzeln 
wahrend der ganzen Wachstumszeit konstant halten kann. 
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In der vorliegenden Arbeit wird nur die Situation erdrtert, die 
sich bei Variierung der Nahrstoffkonzentration u ergibt. Wie die 
Verhaltnisse liegen, wenn nicht wu, sondern 2 variiert wird, soll 
erst in spateren Arbeiten behandelt werden. 

Die Situation bei Mikroorganismen. Die folgenden Betrachtungen 
beschranken sich zunachst auf Bakterien, da hier die Zusammen- 
hange am leichtesten zu durchschauen sind. 

Bezeichnet man die Keimzahl bzw. den Ertrag an Organismen- 
masse mit y und die Zeit mit ¢, so wird die Bakterienvermehrung in 
Abhangigkeit von der Nahrstoffkonzentration u durch die Differential- 
gleichung 


— = yw(u) (1) 


beschrieben (Kaltofen, 1960). An und fir sich kann wu in dieser 
Gleichung eine beliebige Funktion der Zeit sein. Wie bereits 
angedeutet, soll hier jedoch nur der Fall ins Auge gefasst werden, 
dass wu im ganzen betrachteten Zeitintervall konstant bleibt. 
Selbstverstandlich ist Gleichung (1) nur dann giltig, wenn sich alle 
ibrigen Umweltbedingungen mit Fortschreiten der Zeit nicht andern. 

Ist uw = 0, so ist auch wW = 0; denn bei vollstandigem Fehlen eines 
lebenswichtigen N&ahrstoffes vermehren sich die Organismen nicht 
mehr. Mit wachsender N&hrstoffkonzentration wird dann die Bak- 
terienvermehrung zundchst geférdert und danach wieder gehemmt. 
Entsprechend steigt die Funktion w(w) mit wachsendem wu zun&chst 
an, um danach wieder abzufallen. 

Es erhebt sich nun die Frage, wie sich der Ertrag y als Funktion 
der Nahrstoffkonzentration wu verhalt, wenn (1) gilt. 

Liegt zur Zeit ¢ = 0 der Ertrag yo vor, so lautet die Lésung der 
Gleichung (1) 


y= Yo eve. (2) 
Durch Differentiation der Funktion (2) nach u ergibt sich 
oY = yt y’ (uje™™. 3) 
Ou 


Wie (3) zeigt, ist dy/du an den Zeitpunkten ¢ > 0 positiv, gleich 
Null oder negativ, je nachdem ob dy/du = w’(u) positiv, gleich Null 


oder negativ ist. 
Zu sehr interessanten Ergebnissen gelangt man, wenn man (3) 


nochmals nach wu differenziert. Man erhalt so 
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a = yo tlt (’(u))? +o (u)] eh, (4) 
wobei natiirlich w’’(u) =d2y/du? ist. Im folgenden sei nur der 
aufsteigende Ast der Ertragskurve betrachtet, d.h. der Bereich der 
Variablen u, in dem dy/du positiv ist. Wie soeben gezeigt, ist 
letzteres da und nur da der Fall, wo dw/du positiv ist. 

Aus (4) entnimmt man, dass 0?y/du? an den Zeitpunkten ¢ > 0 
positiv, gleich Null oder negativ ist, je nachdem ob der Ausdruck 


A = t(p’(u))? + wu) (5) 


positiv, gleich Null oder negativ ist. Hierbei ist ¢(W’(u))? natirlich 
stets positiv. Die Funktion w’’(w) hingegen kénnte entweder posi- 
tiv oder gleich Null oder negativ sein; diese Frage lasst sich nicht 
ohne weiteres entscheiden. Ist w’’(u) <0, so ist A an der Stelle u 
fir alle hinreichend kleinen ¢ negativ. Jedoch bemerkt man, dass A 
auch bei negativem w’’ (uw) positiv wird, sobald ¢ einen bestimmten 
Wert iiberschreitet. Ist wW’’(w) 2 0, so ist A selbstverstandlich bei 
jedem ¢ > 0 positiv. 

Durch Ausdehnung dieser Betrachtungen auf ganze Intervalle 
ergibt sich der Beweis fiir den folgenden Satz: 

Gilt Gleichung (1), so ist d%y/du? in jedem abgeschlossenen 
Intervall der Variablen u, in dem ws’ (u) > 0 ist und Ww’ (u) und Wy” (u) 
beschrénkt sind, fiir alle t oberhalb eines wohlbestimmten endlichen 
Wertes t, 20 positiv. 

Die Ableitungen wW’({u) und w’’(u) werden wohl stets fir alle 
u 20 beschrankt sein. Insbesondere wird W’(u) fir w—+0 wohl in 
keinem Falle unendlich gross werden. 

Zwecks weiterer Erlauterung der behandelten Zusammenhange sei 
ein Beispiel angefihrt. 

Setzt man in Gleichung (2) yo = 1 und w(u) = 2u -—u?, so erhalt 
man fiir den Zeitpunkt ¢=10 die im oberen Teil der Abb. 1 
dargestellte Ertragskurve y = y(u). Im unteren Teil der Abb. 1 ist 
der Verlauf der Funktion Y = 2u — wu? veranschaulicht. Ersichtlich 
ist d*y/du? durchweg negativ. Trotzdem ist aber 02y/du? in einem 
grossen Bereich der Variablen wu positiv. Erst bei sehr hoher Nahr- 
stoffkonzentration wu wird d®y/du? negativ. Der Wendepunkt der 
Ertragskurve, der bei der Konzentration u = 0.776 liegt, ist in Abb. 
1 durch eine senkrechte Linie gekennzeichnet. 

Zur Zeit ¢= 1 wiirde der Wendepunkt der Ertragskurve bei der 
N&ahrstoffkonzentration u = 0.293 und zur Zeit ¢ = 100 bei u = 0.929 
liegen. Er riickt im Beispiel mit Fortschreiten der Zeit zu immer 
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Abbildung 1. 


hdheren Konzentrationen u. Anders ausgedriickt, der Bereich, in 
dem 07y/du? negativ ist, wird mit Fortschreiten der Zeit immer 
kleiner. Ist uw; <1 beliebig vorgegeben, so ist d?y/du? fir alle 
t > $(1 -—u;)? im ganzen Bereich 0 < uv S w; positiv. Man erkennt 
dies sofort, wenn man die Ableitungen w’(u) = 2 — 2u und w”’ (u) = 
— 2 in Gleichung (4) bzw. (5) einsetzt. 

Die vorstehenden Ausfiihrungen, aus denen unter anderem hervor- 
geht, dass d*y/du? positiv sein kann, beziehen sich zundchst nur 
auf Bakterien. Jedoch diirften diese Ausfiihrungen dem Sinne nach, 
teilweise vielleicht mit gewissen EKinschrankungen, auch fiir Pilze, 
Algen und andere Mikroorganismen gelten. Dariiber hinaus wird 
0?y/du? oftmals auch dann positiv sein, wenn die Tochterindividuen 
nicht—wie im vorangehenden vorausgesetzt—standig genau die 
gleichen Umweltbedingungen vorfinden wie die Ausgangsindividuen, 
wenn sich jedoch die jeweils vorliegende Nahrstoffkonzentration u 
mit Fortschreiten der Zeit nicht Andert. 

Zur Erganzung noch einige Hinweise fiir kiinftige experimentelle 
Arbeiten auf dem betrachteten Gebiet. 

Es ist selbstverstandlich méglich, den oben formulierten Satz ex- 
perimentell direkt nachzupriifen. Jedoch erscheint eine solche 
Nachpriifung wenig sinnvoll. Findet man bei der Versuchsdurchfih- 
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rung in jenem Bereich der Variablen u, in dem dy/du positiv ist, 
keine positiven Werte fiir 0?y/du?, so kann dies nur daran liegen, 
dass die Versuchsdauer zu kurz bemessen war oder dass irgend- 
welche Fehler unterlaufen sind, z. B. dass die jeweils vorgesehenen 
Nahrstoffkonzentrationen wu nicht standig eingehalten wurden. 

Viel wichtiger als eine direkte experimentelle Nachpriifung des 
angefiihrten Satzes ist die experimentelle Bestimmung der Funk- 
tion w =w(u). Kennt man den Verlauf dieser Funktion, so kann 
man den Ertrag y und die Ableitungen dy/du und 0?y/du” fiir jede 
Konzentration wu und fir jeden Zeitpunkt ¢ berechnen. Die Funk- 
tion & ist gewissermassen ein Bestandteil der Reaktionsnorm des 
jeweiligen Organismus. 

Aus den Versuchsdaten erhalt man w nach der Formel 


iefuye log y(u, : — log yo (6) 
Fir y(u, ¢) ist in (6) der zur Zeit ¢ bei der Nahrstoffkonzentration u 
ermittelte Ertrag einzusetzen. 

Voraussetzung dafiir, dass w(u) nach (6) richtig bestimmt wird, 
ist natirlich, dass die jeweiligen Umweltbedingungen bis zum 
Zeitpunkt ¢ konstant gehalten werden. Zur Kontrolle kann man die 
Bestimmung an zwei verschiedenen Zeitpunkten ¢;, ¢2 durchfiihren. 
Zu beiden Zeitpunkten miissen sich die gleichen Resultate ergeben, 
d.h. es muss 


log y(u,t1)-log yo _ log y(u, t2) — log yo 7 
ty ~ to ( ) 


sein. Ist ¢; < ¢g und ist 


log y(u, t1) -— log yo e log y(u, te) — log yo 
ty to ; 


so besteht der Verdacht, dass die Nahrstoffkonzentration im Laufe 
der Versuchszeit abgesunken ist oder dass sich wachstumshemmende 
Stoffwechselprodukte im Kultursubstrat anreicherten. Ist an- 
dererseits 


log y(u, 1) — log yo < log y(u, ta) — log yo 


fi E ; 


80 kann dies ebenfalls auf einem Absinken der Nahrstoffkonzentra- 
tion beruhen. Liegt namlich anfangs eine sehr hohe Nahrstoffkon- 
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zentration vor, so kann ein Absinken der Konzentration durchaus 
vermehrungsférdernd wirken. 

Wenn das Bestehen der Gleichung (7) auch noch keine absolute 
Garantie fir die korrekte Bestimmung der Funktion Ww (u) ist, so ist 
doch bei Bestehen der Beziehung (7) mit einer gewissen Wahr- 
scheinlichkeit anzunehmen, dass w(u) richtig bestimmt wird. 

Bringt man eine geringe Menge Bakterien in ein keimfreies 
flissiges N&hrsubstrat, dem im Laufe der Versuchszeit weder 
Stoffe zugefiihrt noch entnommen werden, so diirften die durch 
fortgesetzte Teilung entstehenden Individuen eine gewisse Zeitlang 
praktisch die gleichen Umweltbedingungen vorfinden wie die Aus- 
gangsindividuen. Sp&terhin, nach starker Vermehrung der Organis- 
men, wird dies jedoch nicht mehr der Fall sein. Bei der Bildung 
von Bakteriensubstanz werden ja Nahrstoffe verbraucht. Ausserdem 
werden durch gréssere Bakterienpopulationen auch entsprechend 
grosse Mengen von Stoffwechselprodukten ausgeschieden. Will man 
also die Funktion W(w) mit Hilfe ‘‘stehender’’ Nahrlésungen be- 
stimmen, so ist die Bestimmung tunlichst vorzunehmen, ehe sich die 
Organismen zu stark vermehrt haben. 

Vielleicht eignet sich auch der von A. Novick und L. Szilard 
(1950) beschriebene Chemostat oder ein ahnliches Gerat zur Be- 
stimmung der Funktion w(u). Beziglich der Theorie solcher Apparate 
sei auf die Ausfiihrungen C. C. Spicers (1955) verwiesen. 

Die Situation bei héheren Pflanzen. Ebenso wie das Wachstum 
von Mikroorganismenpopulationen beruht auch das Wachstum der 
hdheren Pflanzen auf Zellteilung bzw. auf Selbstverdopplung or- 
ganischer Strukturen. Dies hat zur Folge, dass der Ertrag sowohl 
bei Mikroorganismen als auch bei héheren Pflanzen unter bestimm- 
ten Voraussetzungen mit Fortschreiten der Zeit exponentiell bzw. 
ann&hernd exponentiell zunimmt. 

Indessen sind bei hdheren Pflanzen noch gewisse Besonder- 
heiten zu beachten. So haben viele hdhere Pflanzen die Eigen- 
schaft, im Laufe der Zeit zu reifen und dabei mehr oder weniger 
vollstandig abzusterben. In dem Zeitabschnitt, in dem Reife- 
und Absterbevorgange starker in Erscheinung treten, kann natiirlich 
nicht davon die Rede sein, dass der Ertrag mit Fortschreiten der 
Zeit exponentiell wachst. 

Allerdings sind Reife- und Absterbeerscheinungen meist erst 
nach Ablauf einer gewissen Zeit in starkerem Masse zu beobachten. 
Im Anfangsteil der Wachstumszeit sind die Reife- und Absterbe- 
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vorginge meist so schwach ausgepragt, dass man sie unbedenklich 
vernachlassigen kann. Dabei kann der Zeitabschnitt, in dem dies 
unbedenklich ist, ziemlich lang sein. Viele Pflanzenarten, zum 
Beispiel viele Graser, treten im ersten Jahr nicht in die generative 
Phase ein, sondern sie verharren mehr oder weniger in der vegeta- 
tiven Phase. Bei solchen Pflanzenarten kann sich der Zeitab- 
schnitt, in dem der Ertrag mit Fortschreiten der Zeit annahernd 
exponentiell wachst, tiber mehrere Monate erstrecken. 

Ahnlich wie bei Mikroorganismen kann man auch bei hdéheren 
Pflanzen in jenem Zeitabschnitt, in dem der Ertrag y mit Fort- 
schreiten der Zeit ¢ exponentiell zunimmt, das Wachstum dy/d¢ 
dem Ertrag proportional setzen. Dariiber hinaus hangt dy/dé bei 
hdheren Pflanzen wie bei Mikroorganismen noch von der N&hrstoff- 
versorgung ab. Jedoch geht es bei héheren Pflanzen nicht an, dy/ 
dt einfach als Funktion der Nahrstoffkonzentration u zu betrachten. 
Bei Mikroorganismen erscheint dies insofern berechtigt, als bei 
diesen jede einzelne Zelle unmittelbar der Aussenkonzentration u 
ausgesetzt ist. Die weitaus meisten Zellen der héheren Pflanze 
dagegen kommen mit der Nahrlésung nicht direkt in Beriihrung. Sie 
sind nur der ‘‘inneren’’ Nahrstoffkonzentration ausgesetzt. Dem- 
gemass erscheint es verniinftig, das Wachstum dy/d¢ bei héheren 
Pflanzen als Funktion der ‘‘inneren’’ Nahrstoffkonzentration bzw. 
des relativen Nahrstoffgehaltes der Pflanze zu betrachten. Der 
relative N&hrstoffgehalt werde mit v bezeichnet. Die Grésse » 
ihrerseits ist natiirlich irgendwie von der fusseren Nahrstoffkon- 
zentration u abhangig. Darauf wird im folgenden noch eingegangen. 

Es kénnte der Einwand erhoben werden, dy/dt hinge auch bei 
Mikroorganismen nicht unmittelbar von u, sondern von v ab, Streng 
genommen mag dies tatsachlich der Fall sein. Jedoch begeht man 
keinen grossen Fehler, wenn man unterstellt, dass dy/dt bei Bak- 
terienpopulationen usw. direkt von u abhangt. Bei hdheren Pflan- 
zen dagegen fiihrt die Annahme, dass dy/dt unmittelbar von der 
momentanen Aussenkonzentration uw abhangt, zu krassen Wider- 
spriichen mit Erfahrungstatsachen (Kaltofen, 1956). 

Auf Grund der vorstehenden Erérterungen kann man bei hdheren 
Pflanzen in dem Zeitabschnitt, in dem Reife- und Absterbevorgange 
noch nicht starker in Erscheinung treten, 


Ss =y ¥(e) (8) 


setzen (Kaltofen, 1957b, 1958). 
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a 


In welcher Weise hangt nun v von u ab ? Bezeichnet man den 
absoluten Nahrstoffgehalt der Pflanze mit N, so ist (2) = N(t)/y(t). 
Und setzt man dN/dt = N’(t) und N(0) = No, so ist ferner 


N(t) = No + [ N’(t)d¢. 
0 


Die Grésse N’(t), die Geschwindigkeit der Nahrstoffaufnahme, 
hangt nun unmittelbar von der momentanen Niahrstoffkonzentra- 
tion uw ab. Diese Abh&ngigkeit sei durch die Funktion 9(u) aus- 
gedriickt. Ausserdem diirfte N’(¢) der resorbierenden Wurzelober- 
flache W proportional sein. Die Wurzeloberflache ist natiirlich vom 
Ertrag nicht unabhangig. Damit ergibt sich ftir v der Ausdruck 

Cae ag 
y 


al t 
-*, | W(y)p(u)de. (9) 


In diesem Ausdruck kann w an und fiir sich eine beliebige Funktion 
der Zeit sein (Kaltofen 1957b, 1958). Hier wird jedoch nur der Fall 
betrachtet, dass sich die 4ussere N&hrstoffkonzentration wu mit 
Fortschreiten der Zeit nicht 4ndert. 

Setzt man den Ausdruck (9) in (8) ein, so erhalt man die Integro- 
Differentialgleichung 


Y yy (4s + , f Wow) (10) 


Da (u) nicht von ¢ abhangt, wenn sich w mit Fortschreiten der 
Zeit nicht andert, kann man 9(w) auch vor das Integral schreiben. 

Die Gleichung (10) allgemein zu lésen diirfte auf grosse Schwierig- 
keiten stossen. Fiir gewisse Spezialfalle kann man jedoch Lésungen 
angeben. Setzt man etwa W(v) = av und W(y) = by, wobei a und 6 
geeignet gewahlte positive Konstanten sind, so ergibt sich aus 
(10) die Gleichung 


dy ; 

a= aNo + abp(u) i ydt, (11) 
0 

die leicht zu lésen ist. Selbstversténdlich entspricht aber die 

Funktion Y = av nur dann einigermassen der Wirklichkeit, wenn v 

nicht zu hohe Werte annimmt. Bei extrem hohem v ist d¥/dv sicher 

negativ. 
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Als Lésung der Gleichung (11) findet man, wenn zur Zeit ¢ = 0 
der Ertrag 7 vorliegt und ~(u) > 0 ist, die Funktion 


a yo + va No eVabpuyt , 1 Yo _VaNo go N SORE ay 
2 VbH(u) 2 Vv 5p(u) 


Um zu sehen, ob 0?y/du? positiv, gleich Null oder negativ ist, 
kénnte man nun wieder daran denken, die Funktion y zweimal nach 
u zu differenzieren. Der Ausdruck, der sich auf diese Weise fir 
d7y/du? ergibt, lasst sich jedoch nicht so leicht diskutieren wie 
der entsprechende Ausdruck bei den Bakterien. Es sei deshalb nur 
ein Beispiel durchgerechnet. 

Ist @=1, 651, No= 1, yo = 1, ¢= 10cund ou) = log( +1) s6 
erhaélt man nach (12) die in Abb. 2 dargestellte Ertragskurve y = 
y(u). Wie die Abbildung zeigt, ist 0?y/du? im vorliegenden Beispiel 
positiv. Und dies ist der Fall, obwohl d?¥/dv? identisch gleich 
Null und d?~/du? negativ ist. Aus Stetigkeitsgriinden darf man 
annehmen, dass 0*y/du? unter gewissen Voraussetzungen auch 
dann noch positiv bleibt, wenn d?¥/dv? negativ wird. 

Natiirlich besagt das angefiihrte Beispiel nicht, dass 0?y/du? 
bei héheren Pflanzen stets positiv ist. Immerhin ist aber auf Grund 
des Beispiels zu vermuten, dass 07y/du? bei héheren Pflanzen 
wenigstens in bestimmten Fallen positive Werte annimmt. 

Kine Bestatigung findet diese Vermutung durch Wasserkultur- 
versuche F. W. Parkers (1927) mit Roggen. In diesen Versuchen 
wurde durch Steigerung der Phosphatkonzentration von 0.05 auf 
0.10 mg PO," pro Liter teilweise ein héherer Ertragszuwachs 
erzielt als durch Steigerung der Konzentration von 0.00 auf 0.05 
mg PO," pro Liter. Dies bedeutet, dass d?y/du? teilweise positiv 
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Abbildung 2. 
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war. Kine weitere Bestitigung lieferte ein eigener Versuch mit 
Spinat. In diesem hatte die Erhdéhung der Nahrlésungskonzentra- 
tion von 0.2%. auf 0.4%. einen grésseren Ertragszuwachs zur Folge 
als die Erhdéhung der Konzentration von 0.0%} auf 0.2%. 07y/du" 
war also auch hier positiv. Die Durchfiihrung weiterer Versuche 
ware allerdings sehr erwiinscht. Es gilt, naher zu erforschen, 
unter welchen Bedingungen 07y/du? positive Werte annimmt. 

Wahrscheinlich ist 0?y/du? bei hdheren Pflanzen nicht von 
Beginn der Wachstumszeit an, sondern erst von einem spateren 
Zeitpunkt an positiv. Dabei diirfte der Wendepunkt der Ertrags- 
kurve mit Fortschreiten der Zeit zu immer hdheren N&hrstoffkon- 
zentrationen wu ricken. Am ehesten treten positive Werte der 
Ableitung 0*y/du? zweifellos im Bereich niedriger Konzentra- 
tionen wu auf, 

Nachdem nunmehr die Zusammenhange in demjenigen Zeitab- 
schnitt einigermassen geklart sind, in dem Reife- und Absterbe- 
vorgange noch nicht starker in Erscheinung treten, seien die 
Betrachtungen auf die ganze Wachstumszeit einjahriger Pflan- 
zenarten ausgedehnt. 

Dem allm&hlichen Reifen und Absterben einjaéhriger Pflanzen 
kann man dadurch Rechnung tragen, dass man y auf der rechten 
Seite der Gleichung (8) mit einer Funktion A = A(¢) multipliziert 
(Kaltofen, 1958). Man erhalt so die Gleichung 


dy 


ae = YY). (13) 


Hierbei fallt A mit Fortschreiten der Zeit monoton ab, Zur Zeit der 
Totreife wird X=0. Ausserdem kann man (0)=1 setzen. A(z) 
gibt, grob gesprochen, denjenigen Anteil an Pflanzenmasse an, 
der zum Zeitpunkt ¢ noch neue Pflanzensubstanz produziert. Fir 
v ist in (13) wieder der Ausdruck (9).einzusetzen. 

Es zeigt sich nun, dass 0*y/du? auch bei Giltigkeit der Gleichung 
(13) positiv sein kann. Auf die Durchrechnung eines entsprechenden 
Beispiels sei hier verzichtet. Bei friiheren Berechnungen, die auf 
Gleichung (13) basierten, ergab sich unter anderem, dass der 
Ertrag y zur Zeit der Totreife bei der Nahrstoffkonzentration wu = 2 
rund dreimal so hoch lag wie bei der Konzentration wu = 1 (Kaltofen, 
1958). Dies bedeutet, dass 07y/du” positiv war. Allerdings war 
dies nicht vom Zeitpunkt ¢ = 0 an, sondern erst von einem spateren 


Zeitpunkt an der Fall. 
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Es ist also durchaus nicht so, dass 07y/du” nur dann positive 
Werte annehmen kann, wenn der Ertrag mit Fortschreiten der Zeit 
exponentiell oder annahernd exponentiell wachst. Bedeutungsvoll 
fir das Auftreten positiver Werte der Ableitung d*y/du? ist jedoch 
die Tatsache, dass das Wachstum dy/dt um so groésser ist, je 
grosser y zum Zeitpunkt ¢ ist. Diese Tatsache kommt sowohl in 
Gleichung (8) als auch in Gleichung (13) zum Ausdruck. Da die 
gleiche Tatsache unter Umstaénden auch ein exponentielles bzw. 
annahernd exponentielles Wachstum zur Folge hat, ist d7y/du? oft 
dann positiv, wenn der Ertrag mit Fortschreiten der Z eit exponentiell 
bzw. annahernd exponentiell ansteigt. 


Schlussbemerkungen. Aus den vorstehenden Untersuchungen 
geht unter anderem hervor, dass 0*y/du? im aufsteigenden Ast der 
Ertragskurve nicht negativ zu sein braucht, sondern durchaus auch 
positiv sein kann. Dieses Ergebnis steht im Widerspruch zu der 
Ansicht, der ‘‘Ertragszuwachs’’ dy/du falle im aufsteigenden Ast 
der Ertragskurve mit wachsendem wu stets ab. Ihren Ausdruck findet 
diese Ansicht namentlich im sogenannten ‘‘Gesetz vom abnehmenden 
Ertragszuwachs’’, das oft als ein allgemeingiiltiges, die ganze 
organische Welt beherrschendes Naturgesetz hinges tellt wird 
(Mitscherlich, 1952; Scheffer, 1952, S. 397; Lundegardh, 1957; 
Schmalfuss, 1958, S. 228). 

Die Ansicht, dass der ‘‘Ertragszuwachs’’ im aufsteigenden Ast 
der Ertragskurve mit steigender Nahrstoffversorgung stets abnehme, 
hat sich auf Grund gewisser experimenteller Befunde herausgebildet. 
Hierbei handelt es sich aber fast ausschliesslich um Versuche, in 
denen die vor oder bei Versuchsbeginn verabreichte Nahrstoffgabe 
z variiert wurde. Derartige Versuche kénnen wohl zeigen, wie sich 
der “‘Ertragszuwachs”’ dy/de bei wachsendem z verhialt, sie kénnen 
jedoch kaum etwas dartiber aussagen, wie sich der ‘‘Ertragszu- 
wachs”’ dy/du bei steigendem wu verhalt. Ungeachtet dessen hat 
man angenommen, dass dy/du bei steigendem u dasselbe oder 
zumindest ein ahnliches Verhalten zeigt wie dy/dz bei steigendem 
x (Rippel, 1924; Baule, 1925). Anscheinend ist man sich bisher 
noch nicht gentigend tiber die tiefere Bedeutung des Unterschiedes 
zwischen «a und u klargeworden. 


Ebenso wie man Ergebnisse, die sich auf dy/dx bzw. d?y/dx? 
beziehen, nicht schematisch auf dy/du bzw. d?y/du? tibertragen 
darf, ebenso ist es andererseits natiirlich auch nicht angangig, 


ERTRAGSVERLAUF UND NAHRSTOFFKONZENTRATION 209 


Ergebnisse, die dy/du bzw. d*y/du? betreffen, schematisch auf 
dy/dx bzw. d°y/dx? zu iibertragen. Wenn in den vorstehenden 
Ausfiihrungen nachgewiesen wurde, dass 0*y/du? positiv sein 
kann, so ist damit noch nichts tiber d°y/dx? gesagt. Die Verhali- 
nisse, die bei Variierung der Nahrstoffgabe z vorliegen, sollen erst 
in kiinftigen Arbeiten erértert werden. Zum Teil wird dies darauf 
hinauslaufen, sogenannte logistische Wachstumsmodelle zu un- 
tersuchen. 

Doch wie die Ergebnisse kiinftiger theoretischer und experi- 
menteller Untersuchungen auch ausfallen mégen, es ist und bleibt 
eine feststehende Tatsache, dass der ‘‘Ertragszuwachs’’ mit stei- 
gender N&hrstoffversorgung keineswegs in allen Fallen abnimmt. 
Zumindest bei Variierung der N&hrstoffkonzentration uw nimmt der 
““Ertragszuwachs”’ mit steigender Nahrstoffversorgung oft zu. An- 
gesichts dessen kann die Meinung, das ‘‘Gesetz vom abnehmenden 
Ertragszuwachs’’ sei ein allgemeingiiltiges biologisches Grund- 
gesetz, offensichtlich nicht aufrechterhalten werden. 

Zweifellos konnte sich die Uberzeugung von der Allgemeingiiltig- 
keit des ‘‘Gesetzes vom abnehmenden Ertragszuwachs”’ nicht zuletzt 
deshalb so lange halten, weil man sich bisher zu sehr auf die 
empirisch-experimentelle Untersuchung der einschlégigen Fragen 
beschrankt hat. Nicht selten hielt man eine theoretische Durch- 
dringung der Probleme fir iiberfliissig oder sogar fir abwegig. 
Wohl war bekannt, dass der Ertrag an Organismenmasse von den 
N&hrstoffbedingungen abhangt. Auch wusste man, dass der Ertrag 
unter bestimmten Voraussetzungen mit Fortschreiten der Zeit expo- 
nentiell wachst. Man sah jedoch die Abhangigkeit des Ertrages von 
der Zeit einerseits und von den N&hrstoffbedingungen andererseits 
meist nicht im Zusammenhang. Hatte man sich starker um eine 
Synthese der Einzelerkenntnisse und um eine theoretische Durch- 
dringung der Probleme bemiiht, wire man bestimmt schon friiher 
auf die hier erzielten Ergebnisse gestossen. Im Grunde genommen 
beruhen diese Ergebnisse doch nur darauf, dass zwei wohlbekannte 
biologische Phanomene, die Selbstverdopplung organischer Struk- 
turen und die Umweltabhangigkeit aller Lebensvorginge, in Zusam- 
menhang gebracht wurden. Ubrigens kann man bei Beachtung des 
natiirlichen Zusammenhangs zwischen diesen beiden Phanomenen 
schon durch ganz elementare Betrachtungen zeigen, dass der 
“‘Ertragszuwachs’”’ bei Bakterien mit steigender N&hrstoffkonzen- 
tration u ansteigen kann (Kaltofen, 1960). 
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ERRATA 


In the paper: ‘‘Life, Information Theory, Probability, and Physics” 
by N. Rashevsky, Vol. 22, pp. 351-364, the equation on the third 
line from bottom of page 356 should read: 


fim = logg nj. 


Equation (7) on page 356 should read: 


m 
} logs nj = Mp. 
1 


Equation (9) on page 357 should read: 


3 Gm = ™q- 


t=1 


On the third line following equation (10) on page 357 delete the 
word ‘‘from.”’ 

At the end of the fifth line following equation (11) on page 357 
the expression should be: 


m/logem. 
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ANNOUNCEMENT 


PAN AMERICAN UNION’S REQUEST FOR 
OLD SCIENTIFIC JOURNALS 


Under a new program organized by the Division of Science Devel- 
opment of the Pan American Union, General Secretariat of the Or- 
ganization of American States, thousands of issues of scientific 
and technical journals are being transferred from the overcrowded 
libraries of United States scientists to the understocked shelves of 
the libraries and information centers of Latin America. The trans- 
fer is being accomplished through the U.S. Book Exchange, Inc., 
in Washington. 

Scientists of the U.S.A. are being asked to donate to the U.S.B.E. 
files of scientific and technical journals which they no longer wish 
to keep in their libraries. The value of the gift plus the cost of 
making it (mailing charges in this case) are income tax deductible. 
Latin American libraries are being informed of the services of the 
U.S.B.E. and urged to join and make full use of the opportunity. 

Further information about the program may be obtained from the 
U.S. Book Exchange, 3335 V. St., N.E., Washington 18, D.C., or 
from the Division of Science Development, Pan American Union, 
Washington 6, D.C. 
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ANNOUNCEMENT 


GRADUATE TRAINEESHIPS IN MATHEMATICAL BIOLOGY 
AND BIOMETRY 


Training programs designed to prepare students in the application 
of statistical and mathematical methods to biological problems, 
particularly those related to health and medical sciences, now ex- 
ist in more than 20 universities throughout the country. Supported 
by training grant funds from the Public Health Service, NIH, these 
programs provide unusual opportunities for careers in teaching, re- 
search, and consultation. Employment opportunities for biome- 
tricians and mathematical biologists are excellent, with the demand 
by governmental and voluntary health agencies, medical research 
and educational institutions, and industry running far in excess of 
the available supply of trained personnel. 

Programs of study are individually designed to lead to doctoral 
degrees, and in special instances, to other academic degrees. 
Traineeship stipends are provided at various levels depending on 
previous education and experience of the trainee and include al- 
lowances for dependents. Substantially full economic support or 
partial support may be provided, depending upon the proportion of 
time spent in training. 

For those unable to train during the academic year, an unusual 
opportunity is provided by a cooperative Graduate Session of Sta- 
tistics in the Health Sciences sponsored by the program directors 
of the participating universities and made possible by a training 
grant from the PHS, NIH. For information concerning available 
stipends and course offerings at elementary, intermediate, or ad- 
vanced levels for the summers of 1961 and 1962, write Dr. Jacob 
E. Bearman, University of Minnesota, Minneapolis, Minnesota. 
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Applicants interested in full time training are encouraged 


write to a 
Miss Emmarie C. Hemphill, Executive Secretary on 
Advisory Committee on Epidemiology and Biometry 
Division of General Medical Sciences a 
> National Institutes of Health = 
Public Health Service sa 


Department of Health, Education, and Eee 7 i ees 
Bethesda 14, Maryland at 


